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Objectives 


In  February  1959,  a  study  was  begun  at  the  National 
Bureau  of  Standards  to  determine  the  thermal  characteristics 
of  an  underground  concrete  shelt-.r  with  simulated  occupancy. 
Experimental  studies  Indicated  that  theoretical  approaches 
were  desirable  in  order  to  provids  methods  for  predicting 
environmental  conditions  in  shelters  of  any  size,  shape, 
occupancy,  thermal  properties,  and  ambient  conditions. 

This  would  reduce  the  need  for  exhaustive  testing  under  an 
empirical  approach.  Tests  in  a  few  representative  shelters 
could  then  be  used  to  verify  the  applicability  of  those 
methods. 

Scone 

The  report  deals  with  only  part  of  an  analytical 
approach,  confined  to  the  mathematical  solution  of  the 
problem  of  transient  heat  conduction  in  shelter  models 
characterized  by  simplified  geometrical  configurations 
and  a  constant  heat  flux  applied  at  the  Interior  boundary 
of  the  heat  conduction  region.  It  does  not  include  the 
parameters  of  ventilation  and  occupancy,  or  the  mechanisms 
of  transfer  from  these  sources  to  the  interior  surfaces. 

This  division  of  the  problem  had  the  important  advantage 
of  isolating  existing  sources  of  error  to  those  lying  in 
the  heat  conduction  region,  and  provided  a  means  of  study¬ 
ing  the  heat  conduction  aspects  separately. 
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Tvpes  of  Solutions 

The  six  machematlcal  solutions  contained  in  the 
report  are  based  on  the  three  models  described  below. 

A.  One-dimensional 

The  shelter  is  separated  into  six -regions,  five  of 
which  (walls  and  floor)  are  individually  considered  to 
be  layers  of  infinite  thickness  and  infinite  surface 
area.  Constant  heat  flux  is  applied  over  the  entire 
surface.  Each  region  is  considered  to  be  eithsr  a 
single  Jtomogeneous  mass  of  earth,  or  a  composite  of 
concrete  and  earth.  The  initial  temperature  may  be 
either  uniform  or  non-uniform.  Four  solutions  are 
based  on  this  model,  and  these  solutions  are  not  appli¬ 
cable  to  the  finite  (roof)  region.  This  region  is  net 
treated  in  the  report. 

B.  Spherical 

The  shelter  is  represented  by  a  spherical  cavity  in 
an  infinite  homogeneous  medium  with  constant  heat  flux 
applied  over  the  surface  of  the  cavity.  The  initial  tem¬ 
perature  is  uniform. 

C.  Cylindrical 

The  shelter  is  represented  by  an  infinitely  long  cylin¬ 
drical  avity  in  an  infinite  homogeneous  medium  with  constant 
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heat  flux  applied  over  the  surface  of  the  cavity.  The 
init5.al  temperature  is  uniform. 

Solutions  based  on  Models  B  and  C  are  applicable  to 
the  entire  shelter  only,  and  not  to  individual  surfaces. 

Interior  surface  temperatures  predicted  by  each 
solution  are  compared  with  those  obtained  experimentally 
in  the  NBS  test  shelter  series. 

The  solutions  in  the  report  were  chosen  in  order  to 
study  the  basic  parameters  listed  below. 

1)  geometry 

2)  composition  of  region 

3)  initial  temperature  distribution 

4)  thermal  properties 

Comparison  with  NBS  test  data  ind-  the  following 

conclusions: 

1)  geometry 

For  small  shelters,  the  cylindrical  model  gives  generally 
good  predictions,  and  provides  a  rough  design  formula  for 
these  cases  in  which  an  estimate  of  the  heat  flux  is 
available.  One-dimensional  models  would  be  mors  suitable 
for  larger  shelters,  and  the  spherical  model  for  deeper 
shelters. 


2)  composition  of  region 

The  Importance  of  the  composite  assumption  depends  upon 
wall  thickness  and  the  degree  of  dissimilarity  of  thermal 
properties  of  walls  and  earth.  For  shelters  of  construc¬ 
tion  similar  to  that  described,  it  should  not  be  neglected. 

3)  Initial  temperature 

Initial  temperature  non-uniformities  in  the  heat  conduction 
region  affect  the  shelter  temperature  appreciably  only  during 
the  fjrst  week.  Therefore  the  initial  temperature  may  be 
treated  as  uniform,  unless  the  first  week  is  of  comparable 
importance  to  the  second,  e.  g.  in  very  hot  climates. 

4)  thermal  properties 

The  values  of  thermal  properties  have  an  important  bearing 
on  the  shelter  temperature .  Curves  showing  the  predictions 
for  several  sets  of  properties  are  shown  in  the  report,  for 
each  solution.  Variations  in  earth  thermal  properties  due 
to  moisture  migration  are  probably  nou  large  under  the  con¬ 
ditions  of  a  typical  shelter  occupancy  period. 
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MATHEMATICAL  ANALYSIS  05*  TEMPERATURE  RISE 
IN  THE  HEAT  CONDUCTION  REGION 
OF  AN  UNDERGROUND  SKELTER 

by  • 

F.  J.  Jerome  Drapeau  and  Nell  V.  Baggette 


ABSTRACT 

Six  analytic  solutions  of  the  heat  conduction  equation 
to  predict  Interior  surface  temperatures  in  an  underground 
concrete  shelter  are  presented  in  this  report.  Four  solu¬ 
tions  employ  one-dimensional,  semi-infinite  geometry  and 
are  applied  to  the  wall  and  floor  heat  conduction  regions 
of  the  small  NBS  test  shelter.  These  solutions  are  not 
applicable  to  the  roof  region.  Two  solutions  use  spherical 
and  cylindrical  geometry,  respectively,  and  are  applied  to 
the  entire  shelter.  These  solutions  are  not  applicable  to 
individual  surfaces.  The  one-dimensional  solutions  were 
used. .primarily  tc  examine  the  importance  of  the  initial 
temperature  distribution,  the  composite  nature  of  the  heat 
c  nduction  regions,  and  in  conjunction  with  the  other  two 
solutions,  the  geometry  of  the  model.  Surface  temperature 
predictions,  for  selected  surfaces,  together  with  the 
observed  performance,  are  presented  for  the  one-dimensional 
solutions  in  tabular  and/or  graphical  form,  covering  a  range 
of  assumed  thermal  properties.  Similar  results  are  pre¬ 
sented  for  the  spherical  and  cylindrical  solutions  in 
graphical  form,  for  two  sets  of  thermal  properties  and  two 
values  of  equivalent  shelter  radius.  The  results  show  that 
the  one-dimensional  solutions  give  high  predictions,  the 
spherical  solution  low  predictions,  and  the  cylindrical 
solution  generally  good  agreement  with  observed  performance 
of  the  NBS  test  shelter. 
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1.  INTRODUCTION 

1.1  Objectives 

In  February  1959 >  a  study  was  begun  to  determine  the 
thermal  characteristics  of  a  family-size  underground  con¬ 
crete  shelter  with  simulated  occupancy. 

The  first  phase  of  this  study  was  experimental,  with 
the  objective  of  obtaining  performance  data  on  the  environ¬ 
mental  factors  of  temperature,  humidity,  ventilation,  and 
neat  exchange  in  an  underground  installation  of  this  type 
during  fourteen-day  periods  of  occupancy  in  summer  and 
winter.  This  objective  was  accomplished  and  the  results 
[1]*  reported  In  March  1961. 

The  schedule  of  test  conditions  from  which  the  NBS 
performance  data  were  obtained  is  presented  in  the  follow¬ 
ing  table.  This  report  will  compare  analytical  predictions 

of  shelter  surface  temperatures  with  experimental  data  from 
NBS  Tests  3,  ^ ,  and  5  only. 

TABLE  1.1.1 

SCHEDULE  OP  NBS  TEST  CONDITIONS  [l] 


Ventilation  Air 


NBS  Test 

Date  at 

Duration 

Internal 

Flow 

Avg.  Dry 

Dew  Point 

Number 

Start 

of  Taut 

Heat  Input 

Rate 

Bulb  Temp. 

ow-dn-yr 

days 

Btu/hr 

CFM 

•F 

•F 

1 

8-13-39 

7 

1)0 

42 

85 

69 

2 

8-28-39 

7 

2500 

0 

— 

— 

3 

9-13-59 

14 

2500 

42 

85 

69 

4 

10-6-59 

14 

2500 

18 

85 

69 

5 

3-25-60 

14 

2500 

18 

35 

33 

1  Figures  in  brackets  indicate  literature  references  at 
the  end  of  this  report. 


The  experimental  phase  indicated  that  theoretical 
approaches  were  desirable  in  order  to  provide  methods  for 
predicting  environmental  conditions  in  a  shelter  of  any 
size  or  shape  under  various  ambient  conditions,  and  thus 
to  reduce  the  need  for  exhaustive  testing  under  an  empiri¬ 
cal  approach.  Tests  in  a  few  representative  shelters 
could  then  be  used  to  verify  the  applicability  of  those 
methods.  One  such  approach  was  to  develop  mathematical 
equations  which  would  predict  air  temperature  and  relative 
humidity  inside  shelters  characterized  by  the  following 
known  parameters: 

a)  initial  earth  temperatures 

b)  thermal  properties  of  earth  and  wall 
materials 

c)  construction  detail 

d)  shape  and  size 

e)  weather  conditions 

f)  ventilation  rate 

g)  internal  heat  sources 

This  objective  would  also  include  a  comparison  of  the 
theoretically  predicted  results  with  the  performance  data 
obtained  in  the  experimental  phase,  and  with  any  other 
data  on  shelters,  when  available. 

When  these  equations  were  adequately  substantiated  by 
experiment  for  sufficient  ranges  of  the  above-listed  para¬ 
meters  .  they  would  be  used  to  develop  a  design  manual, 
including  extrapolated  ranges. 

1.2  Scope 

The  scope  of  this  report  was  confined  to  the  presents 
tion  of  mathematical  solutions  to  the  problem  of  transient 
heat  conduction  in  a  certain  class  of  shelter  models,  and 
a  comparison  of  predicted  shelter  temperatures  with  ob¬ 
served  temperatures  in  the  NBS  test  series  [1].  The  class 
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of  models  considered  were  those  characterized  by  a  single 
space  variable  (linear  distance,  spherical  radius,  or 
cylindrical  radius),  and  a  constant  heat  flux  applied  at 
the  interior  boundary  of  the  heat  conduction  region. 

This  report,  then,  deals  with  only  part  of  the  whole 
problem--the  heat  conduction  region,  and  all  solutions 
begin  at  the  interior  boundary.  They  do  not  include  the 
parameters  of  heat  generation  and  ventilating  air  inside 
the  shelter,  the  mechanisms  of  transfer  from  these  sources 
to  the  surfaces  by  convection,  radiation,  condensation, 
and  evaporation,  or  weather  conditions  above  the  shelter. 
However,  experience  with  prototype  shelters  has  shown  that 
wall  surface  temperature  is  a  good  index  of  effective  air 
temperature,  especially  in  small  shelters,  where  they  are 
about  equal. 

This  division  of  the  problem  was  possible  because  the 
rate  of  hea^flow  to  all  of  the  interior  surfaces  as  a  func¬ 
tion  of  timeraetermined  experimentally  during  each  test  on 
the  shelter  studied  at  the  National  Bureau  of  Standards. 

This  method  had  the  obvious  and  important  advantage  of 
isolating  existing  sources  of  error  to  those  lying  in  the 
heat  conduction  region.  Poor  agreement  between  observed 
and  predicted  results  could  not  be  attributed  to  poor 
assumptions  associated  with  the  parameters  of  heat  genera¬ 
tion  or  the  coefficients  of  heat  transfer  from  air  to 
surface. 

In  summary,  therefore,  this  first  attempt  at  obtaining 
mathematical  solutions  to  an  inherently  complicated  problem 
was  simplified  by  dealing  with  available  data  measured  at 
the  boundary,  instead  of  undetermined  heat  exchange  pro¬ 
cesses  inside  the  shelter.  It  provided  a  means  of  studying 
some  of  the  heat  conduction  aspects  of  this  problem,  such 
as  the  importance  of  a  prescribed  initial  temperature  dis¬ 
tribution  in  the  region  and  the  effects  of  the  thermal  and 
physical  properties  of  the  materials  comprising  the  region. 
It  also  provided  a  means  of  determining  the  applicability 
of  the  linear  heat  flow  theory  (i.e.  temperature  and  posi¬ 
tion  independent  coefficients),  to  a  problem  of  this  kind. 
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1.3  Types  of  Solutions 

Pour  solutions  contained  in  this  report  were  based 
on  the  one-dimensional  model.  In  this  model,  the  heat 
conduction  region  exterior  to  the  wall  and  floor  bound¬ 
aries  of  the  rectangular  parallelepiped  cavity  was  rep¬ 
resented  by  a  layer  of  infinite  thickness  and  infinite 
surface  area.  Each  region  was  considered  to  be  either 
a  single  homogeneous  layer  of  earth  or  a  composite  of 
concrete  and  earth.  The  heat  flux  was  presumed  to  be 
applied  uniformly  over  the  entire  surface  so  that  heat 
flowed  in  only  one  direction,  normal  to  each  plane  sur¬ 
face.  The  finite  (roof)  region  Is  not  treated  in  this 
report. 

An  Infinitely  thick  layer,  such  as  described  above, 
is  called  a  semi-infinite  region.  Thi3  denotes  a  region 
bounded  on  one  side  only,  o  £  x  <«>in  the  homogeneous  case, 
or  -  *  s  x  <  o  and  o  &  x  <«  in  the  composite  case. 

One  solution  was  based  on  the  spherical  model  In  which 
the  heat  conduction  region  was  represented  by  an  infinite, 
homogeneous  mass  of  earth  bounded  Internally  by  a  spherical 
cavity  of  radius  a.  The  heat  flux  was  presumed  to  be 
applied  uniformly  over  the  entire  surface  so  that  heat 
flowed  radially  in  all  directions. 

In  the  cylindrical  solution,  the  heat  conduction 
region  was  also  represented  by  an  infinite,  homogeneous 
earth  medium,  bounded  internally  by  an  infinitely  long 
cylindrical  cavity  of  radios  a.  The  heat'- flux  was  presumed 
to  be  applied  uniformly  over  the  entire  surface  so  that 
heat  flowed  only  radially  to  the  cylinder. 

The  assumption  of  uniform  heat  flux  application  over 
the  entire  surface  of  each  model  permitted  the  mathematical 
treatment  of  the  Cartesian  (x,  y,  z)  model,  the  spherical 
(r,  6,  0)  model,  and  the  cylindrical  (r,  9,  z)  model  as 
simple  models  of  only  one  (x  or  r)  space  variable. 

Each  mathematical  model  represents  a  simplification 
of  the  actual  physical  shelter,  and  the  effects  of  this 
simplification  will  be  discussed  later  in  the  report. 


•t 
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2.  SUMMARY  OP  PROBLEMS 

Six  mathematical  solutions  to  study  temperature 
predictions  in  an  underground  shelter,  as  determined  by 
heat  conduction  into  the  surrounding  wall  and  floor  regions, 
are  presented  in  this  report.  The  boundary  conditions  and 
assumptions  for  each  of  the  heat  transfer  models  are  summar¬ 
ized  for  each  solution.  The  initial  temperature  distribu¬ 
tion  and  the  heat  conduction  model  are  shown  diagrammatically 
for  each  summary. 

SOLUTION  I  —  One-dimensional,  homogeneous,  semi-infinite 
heat  conduction  region,*  uniform  initial  temperature;  constant 
heat  flux  across  the  surface. 


This  solution  was  used  as  a  first  approximation  to  pre¬ 
dict  NBS  shelter  performance. 


SOLUTION  II  —  Spherical  cavity  bounded  by  an  infinite, 
homogeneous  heat  conduction  region;  uniform  initial  tempera¬ 
ture;  constant  heat  flux  across  the  surface. 


This  solution  was  used  as  a  second  approximation  to 
predict  performance. 
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SOLUTION  III  —  Infinitely  long  cylindrical  cavity 
bounded  by  an  infinite*  homogeneous  heat  conduction  region; 
uniform  initial  temperature;  constant  heat  flux  across  the 
surface. 

This  solution  was  used  as  a  third  approximation  to  pre¬ 
dict  performance.  _ 


SOLUTION  IV  —  One-dimensional,  composite,  semi-infinite 
heat  conduction  region;  uniform  initial  temperature; 
constant  heat  flux  across  the  surface. 


This  solution  was  used  to  observe  the  effect  of  the 
composite  region  assumption. 


CO 
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SOLUTION  V  —  One-dimensional,  composite,  semi-infinite 
heat  conduction  region;  prescribed  initial  temperature; 
constant  heat  flux  across  the  surface. 


This  solution  was  derived  on  the  basis  of  a  prescribed 
Initial  temperature  distribution,  and  used  to  observe  the 
effect  of  such  a  distribution. 


SOLUTION  VI  —  One-dimensional,  composite,  semi-infinite 
heat  conduction  region;  prescribed  initial  temperature; 
constant  heat  flux  across  the  surface. 


This  solution  is  similar  to  the  preceding  solution, 
differing  only  in  the  mathematical  expression  for  initial 
temperature  distribution  in  the  earth  conduction  region. 


r 


4 
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3.  EVALUATION  OF  EXPERIMENTAL 
PARAMETER  DATA 

The  purpose  of  this  analysis  was  not  only  to  develop 
mathematical  solutions,  but  also  to  extract  results  from 
those  solutions  for  comparison  with  experimental 
observations.  For  this  reason,  it  was  necessary  to 
establish  numerical  values  for  the  several  parameters 
appearing  in  each  solution.  This  section  deals  with 
til©  selection  of  these  values  for  each  of  the  following 
parameters,  namely 

1)  average  initial  temperature 

2)  properties  of  concrete 

3)  properties  of  earth 

4)  heat  flux. 


3*1  Average  Xnitial  Temperature 

Solutions  I  through  IV  were  baaed  on  the  assumption  that 
the  initial  temperature  distribution  throughout  each  region 
was  constant  (or  sero  for  mathematical  convenience).  There¬ 
fore,  a  single  value  of  temperature  had  to  be  used  for 
computing  numerical  results. 

Since  the  initial  temperature  distribution  for  all 
regions  deviated  to  some  extent  from  a  constant  value  over 
the  domain  in  question  (see  Appendix  A),  the  average  Initial 
temperature  up  was  calculated  for  each  composite  region  on 
the  basis  of  the  following  heat  balance? 

xithC'.(u0-u0l)  =  xspPc2(u03-u0)  (1) 

where  the  subscripts  1  and  2  denote  concrete  and  earth. 

uQ1  and  uoa  were  calculated  by  graphically  evaluating 
the  integrals 


and 


u(x)dx 


u(x)dx 


(2) 


(3) 


using  the  Interface  as  the  reference  point.  The  functions 
u(x)  are  the  Initial  temperature  profiles  shown  in  Appendix  A. 


The  values  of  average  initial  temperature  obtained  by 
this  method  are  listed  in  Table  3*1*1* 

It  will  be  noted  in  the  graphs  for  Solutions  I  and 
IV  that  the  origins  of  the  observed  temperature  curves 
show®  in  Sections  4  through  7  are  not  necessarily  the 
same  as  those  of  the  predicted  curves.  This  is  because 
the  initial  temperature  of  the  observed  curve  is  the 
actual  value  at  the  interior  surface  (see  Appendix  A), 
whereas  the  initial  temperature  of  any  predicted  curve 
is  the  average  value  shown  in  Table  3.1.1. 
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TABLE  3.1.1 

COMPUTED  INITIAL  TEMPERATURES 


Average  Region  Temperature  un,  °F 


NBS 
Test 
Number  ! 

Heat  Conduction 

Region 

North 

West 

South 

East 

Roof 

Floor 

1 

70.8 

68.7 

69.6 

70.1 

72.2 

64.2 

2 

71.2 

69.3 

69.3 

71.0 

72.7 

64.5 

3 

71.7 

69-7 

68.3 

71.0 

67.5 

64.8 

4 

68.7 

67-7 

67.9 

68.0 

66.8 

64.7 

5 

43.9 

45.4 

45.0 

44.0 

42.8 

30.5 

3.2  Properties  of  Concrete 

The  density  of  the  shelter  concrete  was  taken  to  be 
1m-5  lb/ft  .  This  value  was  based  on  the  total  weight  of 
cement  and  aggregates,  approximately  63,200  pounds,  an 
estimated  water  of  hydration  content  of  2,400  pounds,  and 
the  purchased  volume  of  16.75  cubic  yards.  ' 


The  specific  heat  of  mass  concrete  varies  but 
little  with  the  type  of  aggregate  and  richness  of  mix, 
and  for  practical  purposes  is  usually  taken  as  a  constant 
for  a  given  concrete  at  about  0.22  Btu/lb-°F  [2,  p.585], 

A  mean  value  of  thermal  conductivity  was  selected  from 
among  several  references  L3»4,5,6]  as  1.15  Btu/hr-ft-°F, 
based  on  an  estimated  moisture  content  and  a  physical 
composition  fitting  that  of  the  shelter  concrete  described 
above. 

The  thermal  diffusivity  was  determined  from  the 

relation  ax  =  — -l—  ,  the  subscript  (1)  denoting  concrete. 

Pi  ci 


-12- 


and  the  computed  value  was  O.O36  ftVhr.  These  values  of 
thermal  properties  were  thought  to  be  representative  of 
those  in  the  NBS  test  shelter,  and  will  henceforth  be 
referred  to  as  the  "preferred”  concrete  values. 

The  references  consulted  on  this  subject  listed 
widelv  divergent  thermal  property  values  for  concrete  in 
the  140-155  lb/ft3  density  range.  The  values  of  specific 
heat  ranged  from  0.156  to  O.23O,  thermal  conductivities 
from  0.50  to  1.35*  and  calculated  thermal  diffusivitles 
from  0.016  to  0.046.  As  a  result,  minimum  and  maximum 
values  were  also  used  for  analytical  predictions. 

3.3  Properties  of  Earth 

During  the  excavation  process,  representative  samples 
of  earth  were  removed  for  density  and  moisture  content 
determinations.  Three  distinct  earth  strata  were  found 
and  the  average  values  of  density  and  moisture  content 
were  determined  to  be  approximately  100  lb/ft3  and  17 
percent, 

A  value  of  0.17  Btu/lb-°F  was  selected  as  a  repre¬ 
sentative  average  specific  heat  for  the  three  earth  strata 
in  the  dry  state  [7,  p. 33l •  Since  specific  heat  is  known 
to  vary  with  moisture  content,  the  selected  value  was 
corrected  for  17  percent  moisture,  giving  a  value  of  0.29 
Btu/lb-°F  for  specific  heat  of  the  moist  earth  [7,  p.89]. 

A  determination  of  the  thermal  dlffusivity  was  made 
from  experimental  data  obtained  during  the  NBS  test 
series.  The  solution  of  the  heat  conduction  equation  for 
a  one-dimensional,  homogeneous,  semi-infinite  region  with 
a  simple  .periodic  boundary  condition  at  the  surface  yields 
a  convenient  expression  for  determining  the  thermal  diffu- 
sivity,  based  on  the  time  lag  between  the  maxima  of  the 
ground  surface  temperature  wave  and  the  corresponding  wave 
at  any  depth  below  the  surface.  (This  phase-lag  method  is 
more  accurate  than  the  temperature-amplitude  method  because 
the  values  of  amplitude  are  difficult  to  measure  accu¬ 
rately  [8]). 


The  expression  is 

Pxa 

“  =  Wt£‘  » 

where  P  =  period,  hours 

x  =  depth  of  temperature  measurement,  feet 
t  =  time  lag  between  maxima,  hours 

Although  temperatures  were  measured  at  several  depths 
in  undisturbed  earth  near  the  shelter,  only  two  were  use¬ 
ful  for  this  determination  -  those  made  at  the  six-inch 
and  one-foot  depths.  Measurements  made  at  the  two-foot 
depth  exhibited  no  clearly  defined  wave  form, if  any  at 
all.  The  amplitude  of  the  daily  cycle  decreases  rapidly 
with  depth  and  is  almost  completely  damped  out  at  the 
two-foot  depth  [9].  For  this  reason,  thermal 
diffusivity  determinations  in  deeper  ground  are  extremely 
difficult  on  a  diurnal  cycle  basis,  and  must  be  made  on 
averages  of  annual  cycles. 


Ambient  temperatures  were  used  because  ground  surface 
temperatures  were  not  recorded.  However,  ambient  and 
ground  surface  temperature  maxima  occur  in  phase,  for 
practical  purposes.  It  was  observed  that  the  temperature 
maxima  at  the  six-inch  and  one-foot  depths,  during  three 
early  consecutive  days  of  NBS  Test  3,  lagged  the  ambient 
maxima  by  about  ^  1/4  and  8  1/2  hours,  respectively. 

These  phase-lags  each  indicate  a  thermal  diffusivity  of 
0.026  ft3/hr. 

A  similar  determination  was  based  on  data  obtained 
during  three  late  consecutive  days  of  NBS  Test  5,  with 
exactly  the  same  phase-lags  and  results.  No  precipita¬ 
tion  had  occurred  for  at  least  six  weeks  prior  to  the 
start  of  or  during  NBS  Test  3,  whereas  two  feet  of  snow 
had  completely  melted  a  few  days  before  the  start  of  NBS 
Test  5- 
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The  thermal  conductivity  was  determined  from  the 
same  relation  K2  =  Ps  Ca  ctg,  the  subscript  (2)  denoting 
earth,  and  the  computed  value  was  0.75  Btu/hr-ft-°F. 

These  values  of  thermal  properties  are  likely  to 
be  representative  of  at  least  the  upper  two-foot  layer 
of  earth  surrounding  the  NBS  experimental  shelter,  and 
possibly  deeper.  They  will  henceforth  be  referred  to 
as  the  preferred  earth  values.  However,  in  order  to 
account  for  bhe  likelihood  of  other  values  at  greater 

depths,  and  changes  occurring  in  the  upper  layer  during 
the  intervening  times,  minimum  and  maximum  values  were 
again  selected  to  form  a  range  which  would  account  for 
reasonable  variations. 

3.4  Heat  Flux 

Although  the  term  "constant  heat  flux"  is  used 
throughout  the  discussion  of  Solutions  I  to  VI,  inclu¬ 
sive,  the  actual  value  of  the  heat  flux  F0  was  used  for 
each  plotted  value  of  time.  The  heat  flux  used  for  time 
t  was  a  cumulative  average  of  heat  flow  meter  readings, 
averaged  over  the  Interval  from  o  to  t.  This  procedure 
was  intended  to  approximate,  to  some  extent,  the  effect 
of  the  observed  exponential  decay  of  the  heat  flux. 

This  is  discussed  more  fully  in  Appendix  B  which  also 
contains  a  table  listing  the  cumulative  average  values 
used  in  each  solution. 

The  west  and  north  wall  regions  were  selected  for 
analysis  and  graphical  presentation  of  the  results  because 
they  were  considered  to  be  most  and  least  closely  approxi¬ 
mated  by  the  one-dimensional  model.  The  reasons  for  these 
selections  will  be  discussed  in  greater  detail  in  Section 
10.5. 
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4 .  SOLUTION  I 

4.1  Statement  of  Problem 

Semi-infinite  homogeneous  region.  The  heat  flux  F0  is 
constant  across  the  surface  x  =  o.  The  initial  temperature 
is  uniform. 

4.2  Mathematical  Summary 

The  differential  equation  to  be  satisfied  by  the 
temperature  rise  e  is 


90 

at  =  a 

with  the  following 
x  >  o 

x  =  a 

X  — >  CO 


a3e 

•  x  >0  »  *  >0 

initial  and  boundary  conditions: 
t  =  o  9  =  o 

t>o  -  K  J*®  =  F0 

ax 

t  >  o  e  O 


(1) 

(2) 

(3) 

(4) 


The  solution  of  the  Laplace  transformed  equation  is 


9  =  -2—  exp  (-qx) 
qKp 

and  the  inverse  transformation  yields  the  solution  for 
temperature  rise 


(5) 


, 
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TT 
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-  —  er 

2 

2/aP 
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(6) 


For  an  initial  temperature  u0,  the  temperature  at  any 
distance  x  in  the  region  is  given  by 


u  =  u0+  9 


(7) 


At  the  interior  surface  x  =  o,  Eq.  4.2(6)  reduces  to 

at 
TT  J 


0  = 


_  ZF0 


(8) 
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4.3  Assumptions 

This  solution  [10,  p.  753  was  used  as  a  first  approxi¬ 
mation  for  predicting  shelter  performance  for  all  heat  con¬ 
duction  regions  except  the  roof  region. 

The  assumptions  used  in  this  approximation  are  listed 
and  discussed  below. 

1)  geometry  -  one -dimensional  model 

The  NBS  shelter  deviated  from  the  one-dimensional 
model  discussed  in  Section  1.3  in  that  the  shelter 
surfaces  were  finite  in  area  with  the  result  that 
the  solution,  being  one-dimensional,  could  not 
account  for  heat  lost  to  the  corner  regions. 

Except  for  the  floor,  deviation  from  the  one¬ 
dimensional  model  wa3  also  caused  by  the  addi¬ 
tional  boundary  at  the  ground  surface. 

2)  homogeneous  heat  conduction  region 

Since  the  region  was  considered  to  be  homogeneous, 
only  one  thermal  conductivity  and  one  thermal 
diffusivity  appear  in  the  solution,  i.e.  both 
the  concrete  and  earth  media  are  assumed  to  have 
identical  thermal  properties.  These  properties 
are  assumed  to  be  invariant  with  position  and 
changes  in  temperature  and  moisture  content. 

For  more  detail,  see  Section  10.4. 

3)  uniform  initial  temperature 

Figures  1  through  5  of  Appendix  A  show  that  the 
initial  temperature  distribution  in  each  region 
was  not  uniform  for  all  tests.  The  assumed  initial 
condition  4.2(2)  most  nearly  fitted  the  wall  regions 
of  NBS  Tests  1  and  54  but  did  not  fit,  for  example, 
the  wall  regions  of  NBS  Test  3.  In  order  to  ob¬ 
tain  numerical  results,  the  temperature  u  was 
computed  on  the  basis  of  the  average  initial 
temperatures  u0  listed  in  Table  3.1.1. 
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4.4  Results 

Figures  4.4.1  and  .2  show  interior  surface  tempera¬ 
tures  versus  time,  predicted  by  Solution  I.  for  the 
north  and  west  wall  concretion  regions  of  NBS  Test  3, 
together  with  the  observed  performance.  Figures  4.4.3 
through  .6  show  similar  temperature-time  curves  for  NBS 
Tests  4  and  5* 

Curves  1,  2,  and  3  show  predicted  performance  based 
on  an  estimated  range  of  earth  thermal  properties,  in¬ 
cluding  the  preferred  values.  Curve  4  is  based  on  the 
preferred  concrete  properties ,  values  which  would  be 
reasonable  if  the  concrete  predominated  in  the  heat 
storage. 

Table  4.4.1  shows  predicted  surface  temperatures 
for  the  five  semi-infinite  regions  of  NBS  Tests  3 ,  *+> 
and  5,  evaluated  at  elapsed  times  of  one  and  two  weeks, 
together  with  the  observed  temperatures  and  the  resulting 
deviations.  The  predicted  temperatures  were  computed 
using  the  preferred  concrete  properties,  the  values  pro¬ 
ducing  the  best  general  agreement  with  the  observed 
performance . 

4.5  Discussion 

Inspection  of  Figures  4.4.1  through  4.4.6  shows 
that  the  lower  thermal  properties  gave  better  agreement 
for  the  north  wall  region,  and  that  the  higher  properties 
were  better  for  the  west  wall  region.  The  behavior  of 
the  other  three  surfaces,  in  most  cases,  was  similar  to 
that  of  tne  west  wall  region  in  this  respect.  This  can 
be  seen  from  Table  4.4.1  which  shows  that  in  11  out  of 
15  cases,  use  of  the  preferred  concrete  thermal  properties 
produced  agreement  to  within  1.5°F  at  the  end  of  a  two- 
week  period. 
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A  comparison  of  analytical  predictions  with  experi¬ 
mental  data  from  NBS  Test  4  is  included  in  this  section 
to  show  the  effect  on  the  observed  curve  of  an  adverse 
ambient  temperature  change  which  occurred  about  halfway 
through  the  test  and  continued  for  its  duration.  This 
can  be  seen  in  Figures  4.4.3  and  4.4.4  where  the  value 
of  observed  temperature  at  fourteen  days  is  actually 
lower  than  the  value  at  ten  days,  producing  an  observed 
performance  impossible  to  duplicate  with  models  of  this 
kind.  For  this  reason,  it  will  be  discarded  for  dis¬ 
cussion  purposes  in  the  succeeding  solutions. 

Further  inspection  of  the  figures  for  NBS  Tests  3 
and  5  shows  that  the  predicted  curves  have  a  basically 
different  shape  from  the  observed  curve,  which  levels 
off  more  rapidly  with  increasing  time.  The  continuing 
rise  of  the  predicted  curves  is  due  to  the  one-dimensional 
nature  of  the  solution  in  which  the  heat  flow  is  confined 
to  only  one  direction  in  space.  The  shelter  itself  could 
lose  heat  in  three  dimensions,  with  the  additional  heat 
sinks  emanating  from  the  corners  becoming  more  important 
with  increasing  time.  This  resulted  in  an  observed  curve 
which  leveled  off  rapidly,  appearing  to  approach  a  con¬ 
stant  temperature. 

These  observations  suggest  a  spherical  model,  which 
conducts  heat  in  three  dimensions.  This  approximation 
should  accordingly  produce  a  temperature-time  curve  of 
more  reasonable  shape,  since  the  solution  is  bounded  in 
time. 
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5.  SOLUTION  II 


5-1  Statement  of  Problem 

Infinite  homogeneous  medium  bounded  internally  by  a 
sphere  of  radius  a.  The  heat  flux  F0  is  constant  across 
the  surface  r  =  a.  The  initial  temperature  is  uniform. 


5.2  Mathematical  Summary 

The  differential  equation  to  be  satisfied  by  the 
temperature  rise  9  is 


de  _ 
•§t  “a 


_d 

dr 


(**S)] 


r  2;  a 


t  >  o 


(1) 


with  the  following  Initial  and  boundary  conditions: 


r  a  a  , 

•s 

0 

11 

-p 

8  =  0 

(2) 

r  =  a  , 

t  >  0  , 

t 

ll 

0 

(3) 

r  •—»  00 

t  >  0  , 

0-fr  0 

(4) 

The  solution  for  temperature  rise  is 


a2Fn 

e  =  ~tk 


erfc 


r  -  a 
2/at 


exp 


r  -  a 


(5) 


For  an  initial  temperature  u0,  the  temperature  at  any 
distance  r  in  the  region  is  given  by 


u  =  u0  +  9 


(6) 
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At  the  interior  surface  r  ~  a,  5*2(5)  reduces  to 
0  =  {  1  "  exp(  t)  erfc  (t£)}  (7) 


where  T  =  Ta  . 

a 

5*3  Assumptions 

This  solution  is  available  on  p.  248  [10]  and  was 
used  as  a  second  approximation  for  predicting  overall 
shelter  performance* 

The  assumptions  used  in  this  approximation  are 
listed  and  discussed  below. 

1)  geometry  -  spherical  model 

The  NBS  shelter  deviated  from  the  spherical 
model  in  that  part  of  the  system  was  actually 
finite,  bounded  externally  by  the  ground 
surface.  It  also  differed  in  the  obvious 
matter  of  shape,  which  would  alter  the  local 
heat  flow  pattern. 

Since  tiie  solution  depends  only  on  radial  dis¬ 
tance  and  time,  as  discussed  in  Section  1*3, 
it  can  make  no  distinction  among  the  six 
surfaces.  Thus,  it  provides  only  a  single 
temperature  prediction  at  a  given  radius  in 
all  directions. 
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a)  equivalent  spherical  radius 

Since  the  assumed  shape  in  this  model 
is  a  sphere,  it  was  necessary  to 
choose  a  sphere  of  such  radius  as  to 
have  heat  transfer  characteristics 
similar  to  those  of  the  actual  shelter. 
This  was  done  in  two  elementary  ways, 
as  follows: 

l)  by  requiring  the  sphere  to  have 
the  same  surface  area  as  the 
actual  shelter,  thus  maintaining 
the  same  total  heat  loss  rate. 

On  this  basis,  the  equivalent 
radius  is  5.7  feet. 

ii)  by  requiring  the  sphere  to  have 
the  same  volume  to  surface  area 
ratio  as  the  actual  shelter, 
thus  taking  into  account,  to 
some  extent,  the  difference  in 
shape.  On  this  basis,  the 
equivalent  radius  is  4.0  feet. 

2)  homogeneous  heat  conduction  region 

This  assumption  was  discussed  in  Section  4.3. 

3)  uniform  initial  temperature 

Since  this  model  is  spherically  symmetric  in 
every  respect,  it  follows  that  only  one  value 
of  initial  temperature  can  be  used.  There¬ 
fore,  the  six  average  region  temperatures  uQ, 
listed  in  Table  were  arithmetically 

averaged  to  obtain  a  single  overall  initial 
temperature . 


-29- 


5.4  Resulcs 

Figure  5.4.1  compares  the  observed  interior  surface 
temperatures  with  predicted  values  using  Solution  II, 
for  the  two  values  of  equivalent  spherical  radius  des¬ 
cribed  in  the  preceding  section  and  for  two  sets  of  earth 
thermal  properties,  including  the  preferred  values.  Fig¬ 
ure  5*4.2  shows  similar  comparisons  for  NBS  Test  5- 

The  observed  temperatures  shown  in  these  figures  are 
arithmetic  averages  of  the  smoothed  experimental  data  for 
all  surfaces,  including  the  roof. 

The  values  of  heat  flux  F0  used  in  computing  these 
results  were  arithmetic  averages  of  the  cumulative  aver¬ 
ages,  including  the  roof,  listed  in  Table  B.l  of  Appendix 
B. 

5.5  Discussion 

Inspection  of  Figures  5*4.1  and  5-4.2  shows  that 
better  agreement  was  achieved  between  observed  and  pre¬ 
dicted  values  by  using  the  estimated  minimum  values  of 
earth  properties  rather  than  the  preferred  values.  It 
can  also  be  seen  that  the  equivalent  surface  area  radius 
of  5.7  feet  gave  closer  agreement  than  the  volume  to 
surface  area  radius.  For  the  minimum  values,  the  5.7  foot 
radius  gave  agreement  to  within  1°F  or  less  at  t  =  336  hours 
for  each  test,  whereas  the  4.0  foot  radius  provided  agree¬ 
ment  to  about  2°F. 

Of  the  two  choices  for  equivalent  spherical  radius 
described,  the  spherical  cavity  which  provided  ar.  interior 
surface  area  (a  =  5.7  feet)  equal  to  that  of  the  parallel¬ 
epiped  shelter  would  be  expected  to  give  better  agreement 
between  the  predicted  and  observed  surfaces  temperatures 
because  the  test  shelter  had  a  low  aspect  ratio  and  thus 
the  total  amount  of  earth  serving  as  a  heat  sink  was  more 
nearly  comparable  in  the  two  cavities. 

Further  inspection  of  Figures  5-4.1  and  5-4.2  shows 
that  the  shapes  of  the  predicted  curves  for  the  spherical 
model  leveled  off  more  rapidly  than  those  for  Solution  I, 
as  would  be  expected.  The  effect  of  the  ground  surface 
boundary  can  be  interpreted  as  a  reduction  in  the  volume 
of  earth  available  to  the  shelter  for  heat  conduction,  but  may  not, 
in  all  seasons  of  the  year,  represent  a  reduction  in  the  total  heat 
absorption  capacity  ef  the  region. 
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'  SOLUTION  H  -  NBS  TEST  5 
PREDICTED  SURFACE  TEMPERATURES 
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6.  SOLUTION  III 
6.1  Statement  of  Problem 

Infinite  h  ^mogeneous  region  bounded  Internally  by  an 
infinitely  long  cylinder  of  radius  a.  The  heat  flux  F0 
is  constant  across  the  surface  r  =  a.  The  Initial  tempera¬ 
ture  is  uniform. 


6.2  Mathematical  Summary 

The  differential  equation  to  be  satisfied  by  the 
temperature  rise  9  is 


J*=a 

at 


dzQ  ,  i  ae 

— -  +  — 

Br3  r 


r  &  a  ,  t  >  o 


with  the  following  initial  and  boundary  conditions: 


r  >  a  s  t  =  o,  9=o 


r  -  a  ,  t  >  o  ,  -K 


_ 


=  Fr 


r  t  >  o  ,  9  -s»  o 

The  solution  of  the  Laplace  transformed  equation  is 


(1) 

(2) 

(3) 

W 


«  =LM2?h 

w  qKpK,(qa)  (5) 

An  asymptotic  series  expansion  of  Eq.  6.2 (5)  yields  the 
solution  for  temperature  rise  for  small  values  of  time 


2FfJ  r  oat  "p  f  „  r  -  a  (3r  +  a)  r  ,  r  -  a 

9  *  k  l  ~  N  Mrfc  iw  -  1  **•  J 1  erfo  zfst 
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For  an  initial  temperature  Uq,  the  temperature  at 
any  distance  r  in  the  region  is  given  by 

u  =  u0  +  0  (7) 

At  the  interior  surface  r  ~  a,  6.2(6)  reduces  to 

e  =  ^aFo  Ti  j  •!„  1  T|  +  4iT  -  ^  +  t,.}(8) 

K  L  JtT  h  K/tt  32  J 

where  T  =  ~r  . 


6 . 3  Convergence 

The  first  two  terms  of  this  series  solution,  Eq. 
6.2(6),  were  available  on  p.  339  [10],  and  were  tried  as 
a  third  approximation  for  predicting  overall  shelter 
performance.  However,  the  contribution  from  the  second 
term  was  found  to  be  considerable  for  moderate  values 
of  the  Fourier  number  T.  suggesting  that  one  or  more 
higher  order  terms  might  contribute  significantly. 

Two  additional  terms  were  obtained  by  expanding  the 
modified  Bessel  functions  of  the  second  kind,  which  appear 
in  Eq.  6.2(5),  into  asymptotic  series,  performing  the 
indicated  division,  and  taking  the  inverse  transform  of 
each  term. 


The  expression  in  Eq.  6.2(8)  converges  for  small 
values  of  T,  but  will  diverge  for  large  values.  The  change 
from  convergence  to  divergence  appears  to  occur  at  approxi¬ 
mately  T  =  1.  For  values  of  T  in  the  neighborhood  of  1, 
convergence  may  be  quite  slow,  so  that  many  terms  of  the 
series  would  be  needed  to  provide  reasonable  accuracy.  In 
this  case,  four  terms  of  the  series  were  sufficient  to  pro¬ 
duce  accuracy  of  the  order  of  hundredths  of  a  degree  Fahr¬ 
enheit. 

A  large  time  solution  p.  339  [10] ,  would  be  required 
to  produce  convergence  for  values  T>1. 


6.4  Assumptions 


The  assumptions  used  in  this  approximation  are 
listed  and  discussed  below. 

1}  geometry  -  cylindrical  model 

The  NBS  shelter  deviated  even  more  from  this 
model  than  it  did  from  the  previous  spherical 
model.  In  addition  to  the  reasons  given  in 
paragraph  (l^  cf  Section  5.3,  the  cylinder  is 
assumed  to  be  infinite  in  length,  as  well  as 
circular  in  shape.  However,  as  stated  in 
Section  5.5,  this  model  was  introduced  for 
purely  mathematical  reasons. 

Like  the  previous  solution,  this  solution 
provides  only  a  single  temperature  at  a  given 
radius,  in  all  directions  normal  to  the  cen¬ 
tral  axis  of  the  cylinder. 

a)  equivalent  cylindrical  radius 

The  cylindrical  radius  chosen  go  approxi¬ 
mate  the  shelter  was  obtained  by  assuming 
that  the  shelter  was  a  cylindrical  cavity 
whose  length  ( (,)  was  equal  to  the  longest 
Internal  dimension  of  the  shelter  and 
whose  radius'  was  such  that  the  curbed  area, 
2iTa  i,,  was  equal  to  the  total  shelter  sur¬ 
face  area.  This  method  yielded  an  equlva*- 
lent  radius  of  6.0  feet.  Only  the  curved 
area  was  used  to  represent,  the  total  sur¬ 
face  area  because  no  heat  flows  axially  in 
the  cylindrical  model, 

2)  homogeneous  heat  conduction  region 

This  assumption  was  discussed  in  Section  4.3. 

3)  uniform  initial  temperature 

The  same  average  initial  temperature  used  in 
the  previous  solution  was  also  used  here,  for 
the  same  reasons  of  symmetry. 


-35- 


6.5  Results 

Figures  6.5.1  and  6.5.2  show  Interior  surface  tem¬ 
peratures  versus  time,  predicted  by  Solution  III,  for  an 
equivalent  cylindrical  radius,  a  =  6.C  ft.  These  curves 
together  with  the  observed  performance,  were  plotted  fcr 
NBS  Tests  3  and  5  for  two  sets  of  thermal  properties, 
including  the  preferred  earth  values. 

The  observed  temperatures  shown  on  each  figure  are 
the  same  ones  used  in  Section  5.4,  and  the  values  of 
heat  flux  F0  used  in  the  computations  are  the  same  as 
those  used  in  the  previous  solution. 

6.6  Discussion 

Inspection  of  Figures  6.5.1  and  6.5.2  shows  that 
the  preferred  earth  thermal  properties  gave  better 
agreement  between  predicted  and  observed  temperatures 
for  an  equivalent  cylindrical  radius  of  6.0  feet.  This 
combination  gave  agreement  to  within  less  than  0.5°F 
for  both  tests,  for  elapsed  times  of  both  168  and  336 
hours.  The  shape  and  position  of  the  predicted  curves 
are  an  Improvement  over  those  produced  by  Solutions  I 
and  II. 

At  this  point,  further  improvement  dealing  with 
models  characterized  by  one  space  variable  could  prob¬ 
ably  be  achieved  by  treating  the  composite  case,  or  the 
nonuniform  initial  temperature  case.  These  refinements 
would  introduce  considerable  complexity  into  the  mathe¬ 
matics;  however,  treating  the  one-dimensional  model  will 
avoid  some  of  the  complexity,  and  will  give  some  indica¬ 
tion  of  the  magnitude  of  the  effects.  Since  the  compos¬ 
ite  nature  of  the  heat  conduction  region  is  a  fundamental 
property  applicable  to  all  concrete  shelters,  it  will  be 
treated  in  the  following  section. 


ELAPSED  TIME  -  DAYS 


ELAPSED  TIME  -  DAYS 
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7 .  SOLUTION  IV 
7.1  Statement  of  Problem 

semi-infinite  composite  region.  The  heat  flux F.  is 
constant  across  the  surface  x  =  -l  .  The  initial  tempe 
ture  is  uniform. 

7.0  Mathematical  Summary 

The  differential  equations  to  be  satisfied  by  the 
temperature  rise  9  are 

=  a,  da9i  _£  <  x  <  o  »  t  >  o  (1' 

St  3x3 


S&a  =  a0  ^  .  o<x<oo,  t>o 

ano  $r 

with  the  following  Initial  and  boundary  conditions 


_  £  <  x  <  Of  t—  Of  0J—  0 

o  <  X  <  OO  f  t  —  O  f  =  0 

v  —  —  i  t.  >  o  »  -  Kj  -  F0 

3x 

x  —  O  f  t>Of 


x  =  o 


V  —  K 

t  >  o  »  h  ^  -$r 


X  'SsSt  ,  t  >  o  ,  8s-^o 
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The  solutions  of  the  Laplace  transformed  equations  are 

et  =  F°  r  exp  [-  q,U  +  X)  ]-  8  exp  [-  g,U  -  *)]  | 
q,K.p  t  1  +  8  exp  (-2q ^i)  J 


and  .....  £o  —  /  exp  [  -  g x(f  +  xk)j  j 

q1K,p(a  +  1)  l  1  +  8  exp(-  2q^)  J 

The  inverse  transformations  yield  the  solution  for 
temperature  at  any  point  in  the  concrete  medium 


(10) 


u,  =  uo  +  F0T^y  (-l)V  {  ierfc  -  g  ierfc  i&i&ZlL  \ 

1  n=0  2T$  2t£  > 

and  the  solution  for  temperature  at  any  point  in  the  earth 
medium 


(11) 


u3  =  Uo  +  i-£  (l-e)F0T^  y  (-l)nBn  ierfc  > 

K,  Li  i 

2TS 


Li 

n=o 


(12) 


where  T  =  axt 
i3 


At  the  interior  surface  x  =  -  ,e  ,  Eq.  7.2(11)  becomes 

CO 

=  U0  +  F0T^  £(-l)V  (  ierfc  2j  -  8  ierfc  ~  }  (13) 


n=o 


At  the  concrete-earth  interface  x  =  o  ,  either  7.2(11) 
or  (12)  may  be  used,  giving 


u  =  u0  +  (l-B)F0T*  \  (-l)V  ierfc  ~|} 


n=o 


(14) 


* 


I 


-40- 


»Jt< 


7,3  Convergence 

Since  o  <  ierfc(x)  <;  l//v!  for  all  values  of  x,  the 
convergence  of  (-g)n  is  sufficient  to  insure  convergence 
of  Eqs.  7.2(11)  and  (12).  Therefore  the  series  converges 
whenever  |  8  !  <  1  .  This  condition  is  always  satisfied,  as 
can  be  seen  from  an  examination  of  the  definition: 


8  =  £-Z_i  and 
a  -  1 

whence  -  1  <  8  <  1  for  0  <  a  <  00. 


a  = 


^2  fOg 
*T  V  a  S 


7.4  Assumptions 

Similar  solutions  were  published  in  references  [11] 
and  [12],  using  a  different  origin  for  the  coordinate 
system.  This  solution,  Eqs .  7.2(11)  and  (12),  was  derived 
with  the  origin  located  at  the  interface  for  mathematical 
simplicity. 

The  assumptions  used  in  this  approximation  are  listed 
and  discussed  below. 

1)  geometry  —  one-dimensional  model 

This  assumption  is  the  same  as  for  Solution  I 
and  was  discussed  in  Section  k, 3.  It  will  be 
listed  in  the  following  solutions  without 
further  comment. 

2)  composite  heat  conduction  region 

The  heat  conduction  region  was  a  composite  of 
8  Inches  of  concrete  in  perfect  thermal  contact 
with  an  infinite  earth  medium,  so  that  the 
solution  contains  the  thermal  properties  of 
both  media.  This  assumption  will  be  listed 
in  the  following  solutions  without  further 
comment . 


I 


/ 
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3)  uniform  initial  temperature 

This  assumption  is  the  same  as  for  Solution  I 
and  was  discussed  in  Section  4.3. 


7-5  Results 

Figures  7.5.1  through  .4  show  interior  surface 
temperatures  versus  time,  computed  from  Eq.  7.2(13), 
for  the  north  and  west  wall  conduction  regions  of  KBS 
Tests  3  and  5,  together  with  the  observed  performance. 

Curves  1  and  3  show  predicted  performance  based  on 
the  preferred  concrete  thermal  properties  and  two  sets 
of  earth  properties,  including  the  preferred  values. 

Curves  2  and  4  show  performance  based  on  estimated 
’'maximum"  concrete  properties  and  two  sets  of  earth 
properties,  namely,  estimated  maximum  and  preferred  values. 

Table  7.5.1  shows  predicted  surface  temperatures 
for  the  five  semi-infinite  regions  of  NBS  Tests  3,  4, 
and  5y  evaluated  at  elapsed  times  of  one  and  two  weeks, 
together  with  the  observed  temperatures  and  the  re¬ 
sulting  deviations.  The  computations  were  based  on 
the  preferred  concrete  and  earth  thermal  properties. 

Table  7.5.2  is  similar  to  7*5.1,  except  that  the 
predicted  surface  temperatures  were  based  on  the  maxi¬ 
mum  sets  of  thermal  properties. 

A  numerical  example  for  this  solution  is  shown  in 
Appendix  C. 


ELAPSED  TIME -DAYS 


ELAPSED  TSME  -  DAYS 


ELAPSED  TIME -DAYS 


TABLE  7.5.1 
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7.6  Discussion 

Inspection  of  Figures  7.5*1  and.  .3  shows  chat  the 
preferred  thermal  properties  of  concrete  and  earth 
(Curve  1)  gave  better  agreement  for  the  north  wall 
region.  This  region,  as  previously  pointed  out,  was 
selected  as  the  least  typical  in  the  sense  of  conform¬ 
ing  least  to  the  one-dimensional  model.  The  maximum 
values  of  thermal  properties  (Curve  *0  gave  better 
agreement  for  the  west  wall  region,  as  well  as  for  the 
other  three  regions,  in  most  cases.  These  findings  are 
similar  to  those  pointed  out  by  Solution  I,  and  can  be 
seen  from  Table  7.5*2  which  shews  that  in  10  out  of  15 
cases,  the  maximum  values  produced  agreement  to  within 
1.5°F  at  the  end  of  a  two-week  period.  The  preferred 
sets,  on  the  other  hand  (Table  7*5*1),  produced  agree¬ 
ment  to  within  the  same  amount  In  only  2  out  of  15 
eases.  This  lack  of  agreement  for  the  preferred  sets 
of  values  is  due  to  the  one-dimensional  nature  of  the 
solution  which  has  already  been  seen  to  cause  high 
predictions.  However,  it  should  be  recalled  that  the 
purpose  of  this  solution  is  to  examine  the  importance 
of  the  composite  assumption. 

Curves  2  and  3,  based  on  approximately  equal  frac¬ 
tional  increases  in  concrete  values  on  one  hand  and 
earth  values  on  the  other,  indicate,  when  compared 
with  Curve  1,  that  the  two  media  play  about  equal  roles 
in  conducting  heat  away  from  the  shelter  for  two-week 
periods.  However,  since  the  earth  is  generally  less 
conductive  than  the  concrete,  and  becomes  increasingly 
utilized  for  heat  storage  as  time  increases,  it  would 
ultimately  be  expected  to  control  the  heat  flow. 

This  tendency  can  be  seen  from  further  examination 
of  Curves  1,  2,  and  3  in  Figs.  7*5.2  and  which  shows 
that  increasing  the  concrete  properties  lowers  the  tem¬ 
perature  more  than  Increasing  the  earth  properties  in 
the  first  week,  but  less  in  the  second  week.  However, 
the  lover  temperatures  at  the  end  of  two  w«eks  (Curve  3) 
are  not  significantly  lower  than  Curve  2,  and  the  reason 
that  the  greater  mass  of  earth  involved  is  not  a  dominant 
factor  by  this  time  will  be  discussed  in  Section  10.2. 

This  solution  can  be  further  improved,  at  the  expense 
of  some  increase  in  mathematical  complexity,  by  consider¬ 
ing  the  non-uniform  initial  temperature  distribution. which 
actually  existed.  (See  Appendix  A). 
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8.  SOLUTION  V 
8.1  Statement  of  Problem 

Semi-infinite  composite  region.  The  heat  flux  F0  is 
constant  across  the  surface  x  =  •  (  .  Each  medium  has 
initial  temperature  u  =  f(x). 


8.2  Mathematical  Summary 

The  differential  equations  to  be  satisfied  by  the 
temperature  u  are 


and 


at  ax2 


P  <  X  <  o  t  t>o 


du2  =  aP  ,  o  <  x  <  CO  ,  t>o 

at  ax2 


(i; 


(2) 


with  the  following  initial  and  boundary  conditions: 


-  .e  <  x  <  o  ,  t-o,  u,=ax4-b  (3) 

o  <  x  <"  oo  ,  t  =  o  ,  u?  =  c  +  (b-c)  exp  (-mx)  (4) 

x  =  -  l  ,  t  >  o  ,  -  Kj  =  F0  (5) 

*x 

x  =  o  ,  t  >  o  ,  Uj  =  ua  (6) 

x  =  o  ,  t  >  o  ,  K,  =  Kg  ^9  (7) 

ax  gx 


,  u2  ->  c 


X  — >  ©o  ,  t  >  0 


(8) 


* 


-so- 


The  solutions  of  the  Laplace  transformed  equations 


are 


-  =  J  F°  2  aKt  1  f 
q,KlP  J  i 


exp  E-q^i  +  x)]  -  P  exp  [-x^U  -  x)] 


1  +  P  exp  (-2q1i,) 


} 


{ 


(b  -  c)  KP  (qs  -  m) 


4-  i  ■  £  |  |  exp  [-g,  ( 2jg  +  x)  ]  +  exp  (q^.)  *> 

1  +  8  exp  (-2q,je)  J 


qjK,  (a  +  1)  (p-opm2) 


/  cr(c-b)-a  ,  r  exp  [  ~g,(2i  +  x)]  +  exp  (g,x)  |  ax  + 

^  q,  p  (a  +  1)  J  ^  1  +  B  exp  ( -2qa j*)  J  +  “ 


and 


-  _  r  Fq  +  aK,  -1  r  exp  [  -qt  (1  •<-_  xk)  3-  p  exp  [  -q1  (l  +  xk)] 


q,Klp  J  1  +  3  exp  C-2qTJt) 


r  q,a  (c-b)-a  ->  r  exp  (-q^k)  +  exp  [-qj(2£  +  xk)]-| 
E  qxP  (cr  +  1)J  t  1  +  B  exp  (-2q,i)  ) 


(b  -  c)(q1  Kta  -  mKs)  ,  exp  (-q^xk)  +  exp  [^(2/  +  xk)3 


f  W  -  ~  r 

qjK^ff  +  l)(p  -  a3ma )  l 


1  +  B  exp  (-2q5i) 


1 - 2l_  |  |  exp  (-tux)  -  exp  (-q^xk)  ]• 


p  -  a2m 


(b-c)  exp  (.x^xk)  +  c 


} 


/ 


(9) 


P 


(10) 
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The  inverse  transformations  yield  she  solution  for 
temperature  at  any  point  in  the  concrete  medium 


CD 

-  7*  (F0  +  aK,  )T*  y(-B)n  ierfc  (2ntl>  +  .  Bierfc  1 ' 

K'  n=o  L  ^3  J 

-  sutd-fOT*  f  (-*)»  [ierfc  x/f  +  ierfo  gj  -,*/«  ] 

nro  L  2T*  21*  J 

. (b-c) ii!>  y(„B)n [.rfC +  erf0 a-iVj ] 

2  „7„  L  2t£  2t4  j 


+  (b-c) 


U+P) 


!<-*>"  { 


errcf  Zz1),?*/*  +  in!  T*1  exp  \(  *YT 


*  2(n+i)  p.  *  “  1  +  erfc  [ 

k  k  J 


£m  skI-i 

+  2 n - L+ax+b 

k  k  J 


2n  -  x£ 

- r~~  +  —  T 

zrt  k 


*H(rJ 


and  the  solution  for  temperature  at  any  point  in  the  earth 
medium 


(F0  +  aK^T^f  (-0)n(  1-B)  ierfc 

n=o  zri 


a£(l-e)T^^  (-8)n  ierfc  2(n+1)  j  ^ 1  +  ierfc  20  +  Wf  1 
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-  (b-c)  I  (-P )n  f" erfc  2(n+l)  +  k x/x  +  erfc  2n  +  kx/x  * 

*-  zr^  zr^  - 


n=o 


+  (b-o)  (-»)“{  erfc£  ^T~~  ♦“**]«*  [(r7T 


+  2(  n+1)  iSL  +  tnx 

k 


+  erfc  [ Zn***/l  + 

-  2T^  k 


exp  [C?Tt 


+  2n  +  mx 


}  +  ¥MirT*-i^H(rJTH 


-  erfc 


/m  T^+  kx/x 
k  2tI 


exp^ilH^T  +  tnx 


T+  (b-c)  erfc  +  c 
}  2Te 


(12) 


where  T  = 

/2 


At  the  interior  surface  x  =  -  x  ,  Eq.  8.2(11)  becomes 


u,  =  -  (F0  +  aK, )  T*f  (-p)n 


n=o 


ierfc  — j-  -  3  ierfc 

T?  Ts 


CO  .  CO 

-  2a*(l-(J)T?y  ( -B)n  ierfc  ~±  .  (I4f})(b-c)  Y  (-3)n  erfc  ££ 
^  2Tt  ^  2T? 

n=0  n  =/*> 


n=o 


(b-c)( l+ft)  Y  ( -B)n  erfc  p“|“  +  ~  jexp£  T 

n~o  ^ 


+  ( 2n+l)  — 
k 


-  aX  +  b 


(13) 
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At  the  concrete-earth  Interface  x  =  o,  either  8,2(ll) 
or  (12)  may  be  used,  giving 


u  =  ZL  (F.  +  aK.)T^y  (-8)nd-a)  ierfc 
K,  2T& 


CO 

-  ai(l-0)T^y  (-g)n  Fierfc  +  ierfc  !1_ 
n=o  L  T*  Tf  J 


,  ,  (1+3) 

-  <b"c>  ~ 


+  (b-c) 


( i+S) 


!<-»"[ 

1=0 


n+1 

erfc  ~T~  +  erfc 
T* 


“*  n+1  fm  i 

erfc  ~r~  +  —  Tf 

.  T*  k 


*H(£f 


T  +  2(  n+1) 


+  erfc 


T  +  2n  ~  j  J-  +  b 


8.3  Convergence 

The  first  two  series  in  Eq.  8.2(11)  are  of  the  same  form 
as  that  of  Eq.  7-2(ll),  and  converge  for  the  reasons  given 
in  Section  7.3- 

The  third  series  contains  complementary  error  functions 
which  are  also  bounded 


o  <  erfc(x)  £  1 

for  all  values  of  x.  Therefore  this  series  al3o  converges 
whenever  |  g  |  1. 
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The  remaining  series  converges  for  all  values  of  T, 
whenever  a  certain  condition  among  the  parameters  3 ,  m, 

1 )  a, ,  and  aa  is  satisfied.  This  condition  can  be  obtained 
by  use  of  the  ratio  test  for  convergence.  The  absolute 
value  of  the  ratio  of  the  (n-fl)k*1  term  to  the  n^h  term 
the  series  is 


|en+1 

1 

erfc 

p2(n+l )-xU  ±  im_ 
.  2T^  ^"T 

M 

exp[ 

^  Jim) 

‘t  +  2(11+1)— 1  .  S-jl 

3n 

erfoT^+lM] 

L  2T*  k  J 

exp 

[(? 

)V 

(1) 


This  reduces  to 


rerf0r2ir4£+  +  41 

L  k  TaJ 

erfcr2«^+  “rTi 

L  2T*  k  J 


(2) 


The  value  of  this  expression  must  be  less  than  one,  in  the 
limit  of  large  n,  for  convergence.  The  bracketed  term 
approaches  one  for  large  n,  and  thus  the  ratio  of  terms 
reduces  to 

|p|exp(2Tr)  (3) 

The  condition  for  convergence,  then,  is  that  (3)  be  kept 
less  than  one.  The  smaller  the  value  of  T,  the  faster  is 
the  convergence,  and  for  T  corresponding  to  about  14  days, 
only  a  few  terms  are  required. 

Eq,  8.2(12)  is  similar  to  (11),  with  two  additional 
terms  which  are  not  summed.  These  terms  vanish  for  large  T. 
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8,4  Assumptions 

A  solution  of  this  type  was  not  available  in  any  of 
several  apposite  references  examined,  and  was  derived, 
as  in  the  previous  solution,  with  the  origin  of  the  coordinate 
system  located  at  the  interface. 

The  assumptions  used  in  this  approximation  are  listed 
and  discussed  below. 

1)  geometry  —  one-dimensional  model 

2}  composite  heat  conduction  region 

3)  non-uniform  initial  temperature 

In  the  previous  solutions,  the  Initial 
temperature  distribution  throughout  each 
region,  e.g.,  -  l  <  x^eo,  was  assumed  to 
be  uniform,  and  an  appropriate  average 
value  was  used  to  represent  the  observed 
distribution  shown  in  Appendix  A.  The 
figures  in  this  appendix  show  that  the 
initial  temperature  distribution  was 
not  uniform  in  most  cases,  especially 
at  the  beginning  of  NBS  Test  3.  Even 
though  an  initial  condition  such  as  that 
of  Test  3  would  not  be  likely  to  exist 
in  an  actual  situation,  especially  one 
in  which  the  shelter  had  been  unoccupied 
for  a.  long  period  of  time,  the  kind  of 
unusual  distribution  which  did  exist 
afforded  an  opportunity  to  observe  its 
maximum  effect  and  significance. 

In  this  solution,  the  Initial  temperature 
was  approximated  in  the  concrete  medium 
by  a  linearly  increasing  function, 
u,  =  ax  +  b,  Eq.  8.2(3),  and  in  the  earth 
medium  by  the  exponential  function 
Ug  =  c  +  (b-c)  exp  (-mx),  Eq.  8.2(4). 

See  Fig.  8.4.1  and  Table  8,4.0. 
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TABLE  8A.0 

VALUES  OP  CURVE-PITTING  CONSTANTS 
FOR  SOLUTION  V 


Curve-Fitting  Constants 


Test 

Surface 

a 

b 

0 

m 

3 

North 

2.40 

66.6 

78.0 

0.50 

West 

1.65 

65.2 

74.0 

0.60 

South 

1.00 

61.7 

7$.o 

0.50 

East 

1.00 

65.0 

77.0 

0.55 

Floor 

1.40 

63.8 

66.0 

0.60 

4 

North 

0.30 

66. 1 

73.0 

0.40 

West 

0.00 

66.0 

70.0 

0.40 

South 

0.00 

66.0 

71.0 

0.40 

East 

0.00 

65.8 

72.5 

0.30 

Floor 

-1.00 

64.7 

64.7 

0.00 

5 

North 

0.00 

45.2 

41.5 

0.30 

West 

0.00 

46.0 

43.5 

0.30 

South 

0.00 

46.0 

42.0 

0.27 

East 

0.00 

45.6 

40.5 

0.30 

Floor 

1.00 

48.3 

54.0 

0,33 

I - 
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This  figure  partly  reproduces  Fig.  3  of 
Appendix  A,  showing  the  initial  tempera¬ 
ture  distribution  in  the  north  and  west 
wall  heat  conduction  regions,  selected  as 
representative  regions  for  detailed 
analysis  in  this  report. 


8.5  Results 

Figure  8.5.1  shows  interior  surface  temperatures 
versus  time,  predicted  by  Solution  V,  for  the  north  wall 
conduction  region  of  NBS  Test  3 ,  together  with  the  ob¬ 
served  performance. 

Curves  1,  2  and  4  in  Figure  8.5.1  were  based  on  the 
preferred  earth  thermal  properties  while  varying  the  con¬ 
crete  thermal  properties  to  cover  the  range  of  values  men¬ 
tioned  in  Section  3,2. 

Curves  1,  3  and  5  in  Figure  8.5.1  were  based  on  the 
preferred  concrete  thermal  properties  while  varying  the 
earth  properties  to  cover  an  estimated  range  of  values  that 
might  occur  at  the  NES  test  site. 

Figure  8.5.2  shows  similar  temperature-time  curves  for 
the  west  wall  conduction  region  of  NRS  Test  3. 

Figure  8.5.3  shows  a  temperature-time  curve  for  the 
west  wall  conduction  region  of  NBS  Test  5.  The  figure 
shcv?s  predicted  performance,  together  with  the  observations, 
based  on  the  preferred  thermal  properties  only. 

The  values  of  the  curve-fitting  constants  used  in  the 
computations  are  listed  in  Table  8,4,0. 

The  results  of  this  solution  were  obtained  on  the  IBM 
7090  digital  computer  at  the  rate  of  0,35  minutes  per  run. 
Each  run  consisted  of  computing  temperatures  at  four  posi¬ 
tions  in  the  concrete  medium  and  eight  positions  in  the 
earth  medium,  for  the  eight  values  of  time  shown  in  the 
figures,  a  total  of  9 6  calculations. 


SOLUTION  3ZI-N8S  TEST  3 
PREDICTED  SURFACE  TEMPERATURES 

FOR 

NORTH  WALL  REGION 


ELA 


ELAPSED  TIME  ~  DAYS 
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806  Discussion 

Comparisons  of  Curves  1  and  2  in  both  Figures  8.5.1 
and  8.5.2  show  that  an  increase  of  14  percent  in  the  thermal 
diffusivity  of  the  concrete  decreased  the  inside  surface 
temperature  0.5-2F  or  less  throughout  the  14-day  period. 

On  the  other  hand,  comparisons  of  curves  1  and  4  in  Figures 

8.5.1  and  8.5.2  show  that  a  decrease  of  a  little  over  50 
percent  in  the  thermal  diffusivity  of  the  concrete  frr  i  the 
preferred  value  increased  the  temperature  about  1°F  f  r  the 
north  wall  surface  and  about  2°F  for  the  west  wall  surface. 

Comparisons  of  Curves  3  and  5  in  Figures  8.5.1  and 

8.5.2  show  that  an  increase  in  thermal  diffusivity  of  the 
earth  from  0.02  to  0.03  fts/hr  caused  a  change  of  less  than 
0.5°F  at  any  time  in  the  surface  temperature  of  the  north 
wall  and  a  maximum  change  of  about  1°F  in  the  surface  tem¬ 
perature  of  the  west  wall  after  14  days.  Thus  it  appears 
that  the  thermal  properties  of  the  concrete  in  the  composite 
model  affects  the  surface  temperature  of  the  wall3  more 
significantly  than  the  thermal  properties  of  the  earth. 

Comparison  of  Figure  8.5.3  with  Curve  1  of  Figure  7.5.4 
for  Test  5  shows  much  similarity  between  the  two,  especially 
after  the  first  week,  indicating  that  when  the  initial  tem¬ 
perature  is  reasonably  uniform  (See  Fig.  5  of  Appendix  A), 
the  predictions  of  Solution  V  approach  those  of  Solution  IV. 

It  was  pointed  out  in  Section  7.6  that  Curves  2  and  3 
of,  eeg.,  Figure  7.5.2,  cross  at  about  7  days,  showing  the 
tendency  of  the  earth  region  to  succeed  the  concrete  as  the 
dominant  medium.  In  Solution  V,  however.  Curves  2  and  3  of 
Figures  8.5.1  and  8.5.2  do  not  cross,  even  after  14  days, 
indicating  that  the  concrete  remained  dominant  throughout 
this  period.  The  reason  is  that  in  Solution  V  the  concrete 
was  assumed  to  be  initially  colder  than  the  earth,  as  shown 
in  Figure  8.4.1,  and  therefore  dominated  the  heat  flow  for 
a  longer  time  than  usual.  This  medel  serves  to  illustrate 
i,he  fundamental  interdependence  of  the  composite  nature  of 
the  region  and  the  nonuniform  Initial  temperature,  which 
could  only  be  seen  in  a  solution  containing  both  assumptions. 


/ 
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9  •  SOLUTION  VI 
9.1  Statement  of  Problem 

Semi-infinite  composite  region.  The  heat  flux  F0  is 
constant  across  the  surface  x  =  -  /  .  Each  medium  has  initial 
temperature  u  =  f(x). 


9.2  Mathematical  Summary 

The  differential  equations  to  be  satisfied  by  the  temper¬ 
ature  u  are 


77  =  ai~>  -i<x<0,  t>o  (i) 

dt  3x3 

aus  a3us 

7^3  3“.  o<x<co,  t>c  (2) 

with  the  following  initial  and  boundary  conditions: 


l  <  X  <  0 

t  =  0 

=  ax  +  b 

(3) 

O  <  X  <  00 

t  =  0 

ua 

=  b  +  cx  exp  (-mx) 

(4) 

X 

II 

1 

t  >  0 

- 

11 

0 

(5) 

ax 

X  =  0 

t  >  0 

=  us 

(6) 

X  =  0 

t  >  0 

K, 

=  KS^£ 

(?) 

3x 

ax 

X  _>  00 

t  >  0 

—$>  b 

(8) 

V 
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The  solutions  of  the  Laplace  transformed  equations  are 


u  _  f  Fo  +  *1  f  [  -qt  U+; 

qtKap  '  1  + 


-qv  U+x)]-  P  exp  [  -qt  (l-x)] 1 


8  exp  (-2qtje) 


a  •»  ,  exp  [  -qx  ( Zl  +  x)  ]  +  exp  (q,  x) 


{qTpTSfn  }  { 


1  +  R  exp 


} 


{ 


cK, 


_  |  f  exp  [  -^(21  +  x)]  +  exp  (g^x)  \ 

qTK,  (cr+1)  (p-  aatn8)  1  +  6  exp  (-2 qt/) 


-  f  aagCmKa  (qa-m)  r  exp  (21  +  x)3  +  exp  (q,x)  1  ax  +  b 

*UX  (S+l)  (p  I  1  1  +  p  exp  (-2q1l)  }  -4-  (9) 


and 


u#  =  fF°  +  ^  \  /  C1-^]  0XP  [-g*  U  t  xk) ]  •> 

L  giK,p  ^  1+0  exp  (-aq*/,)  ^ 

r  a  1  f  exP  [-^1  (2jC  +  xk)]  +  exp  (~g,xk)  \ 
*“q,p  ( o+l)  *  *■  1+fi  exp  (~2qxt.)  ' 


{ 


cKs 


\  fexp  pg^  (2 A  +  xk)  ]  +  exp  (-g,xk)  'i 


}{ 


q^Kj  ( o+l)  (p  -  asra8)  J  ^  1+fi  exp  (~2q, £)  •* 


A- 


2^sCmKs  (q8-m) 


l-qTK1  (o+l)(p  -  a8m8)8 


H 


exp  [-g^  (21  +  xk)  '3  +  exp  (-g^xk) 
1+0  exp  (-2qxi) 


} 


2ascm 


}  { e!‘p(-q>’4i)  -  } + 


cxBxp(  -mx)  b 
li  ~l  + 
p  -  aams  p 


(10) 


il 

/ 


*-! 


*1 
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The  inverse  transformations  yield  the  solution  for 
temperature  at  any  point  in  the  concrete  medium 


u,  =  ]  I*  y  (-e)nrierto^+1^/-t  .  BiarfciSsa^l 

K>  A  L  2T4  21*  J 


n=o 


.  atfl-W*  y  ♦  ierf^l 

L>  *-  2T®  2T*  j 


n=o 


c  (HP) 
m  2 


exfC(TjT  +  2<n+1)f +  r]} 


c  (1+g) 
m  2 


« 

)  V  .  _nfP  Jln\zm  im  nKT  r2n-x/i  /m  n 

-L  <-B)  U<tJ  T  -  ^  -  rrra +  tt*J 

n=o 


*4f> +  -  r]} 


*L'V  2Ti. 


J]} 


+  ax  +  b 


(11) 
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and  the  solution  for  temperature  at  any  point  in  the 
earth  medium 


1  n=o 


(2n+l)4kx/jt 
Z?‘ 


.  a,aW  S 

2,  L  ZT*  2Tf  J 


n=o 


-i&pi  + + -M—f^ + f«*] 


n=o 


-<rJT  +  2(n+l)i£  +  mxj  j 

_  £  (Ml  V  (.0)»{[2(^Jt  +  2n£  +  +  T1*] 

n=o 

““[(t/1  +  ^  +  "*]} 

n-o 

*  &?>  •  w  •  fis-ss )'’  •  -i 


-  {»)'  ■  •  tSK©''  •“])•* 


where  T  - 


(12) 


A 

1 

\ 


s 
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At  the  interior  surface  x  =  -  4,  Eq.  9.2(H)  becomes 

Ul  =  X  ("^n[ierfc^l  ~  gicrfc-jJ 

n=o 


-  2ai(l-3)T^  Y  (-0)nierfc — * 

„ S±  2T? 

n~o 


n=o 


S 

5 


i 

i 

% 


i 

4 


It 
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At  the  concrete-earth  interface  x  =  o,  either 
9.2(11)  or  (12)  may  be  used,  giving 


u  =  —OVaKj  1 1  ^  f  (“3)n(l~3)i8rfc^~ 
Ki  L  2T2 


n=o 


,  as 

-  al(L.3)T®  7  (-P)nrierfc£ii  +  ierfcJ-*] 

»«  T*  T*J 


C  (1-fP) 

m  2 


n=o 

p[(f)t  +  «n*l)fi]} 


c  (l; e) 

m  2 


"J  (-8)n{[2(l!)T  +  2n7j=rfof^  +  -FT*] 


^(f>  * 


It 

Ji 


ll 


j 

•  * 


»  J 

i  * 


i  i 


(14) 


4 


V 


1 


h 


\ 
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9*3  Convergence 

Equations  9.2(11)  and  (12)  converge  subject  to  the 
same  condition  discussed  in  Section  8.3,  namely 

|  8  i  exp^ 2“  ^<1  (1) 

9*1*  Assumptions 

This  solution  was  not  available  in  the  references 
examined,  and  was  derived  like  the  previous  one.  It  differs 
from  Solution  V  only  in  the  form  of  the  initial  temperature. 

The  assumptions  are: 

1)  geometry  —  one-dimensional  model 

2)  composite  heat  conduction  region 

3)  non-uniform  initial  temperature 

As  stated  in  Section  8.6,  the  region  beyond  4  feet 
was  expected  to  contribute  to  the  surface  tempera¬ 
tures,  resulting  in  the  need  for  a  better 
approximation  to  the  initial  profile  in  this 
region.  In  this  solution,  the  initial  tempera¬ 
ture  was  approximated  in  the  concrete  medium  by 
a  linear  function  ux  =  ax  +  b,  Eq.  9-2(3), 
as  in  Solution  V,  and  in  the  earth  medium  by  the 
function  u8  =  b  +  cx  exp  (-mx),  Eq.  9.2(4),  shown 
in  Pig.  9*4.1  (See  Table  9.1+.0).  This  figure 
is  completely  analogous  to  Fig.  8.4.1.  It  can  be 
seen  that  the  function  ua  has  a  maximum  correspond¬ 
ing  to  that  suggested  by  the  observed  curve,  and 
subsequently  decays  exponentially  to  the  value  b 
at  infinite  x.  Thi3  is  probably  closer  to  the 
actual  situation  for  a  much  greater  distance  into 
the  earth,  and  should  improve  the  solution  for 
times  greater  than  one  week. 
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TABLE  9.V.0 

VALUES  OP  CURVE- PITTING  CONSTANTS 
POR  SOLUTION  VI 


Curve-Pitting  Constants 


Test 

Surface 

a 

b 

c 

m 

3 

1 .  11  "  1  i 

North 

2.40 

66.6 

5.23 

0.20 

West 

1.65 

65.2 

5.96 

0.30 

South 

1.00 

61.7 

6.85 

0.20 

East 

1.00 

65.0 

5.79 

0.20 

Floor 

1.40 

63.8 

2.99 

0.50 

4 

North 

0.30 

66.1 

2.09 

0.10 

West 

0.00 

66.0 

1.46 

0.15 

South 

0.00 

66.0 

1.57 

C.10 

East 

0.00 

65.8 

1.83 

0.10 

Floor 

-1.00 

64.7 

0.00 

0.00 

5 

North 

0,00 

45.2 

-1.25 

0.15 

West 

0.00 

46.0 

-0.60 

0.15 

South 

0.00 

46.0 

-1.15 

0.15 

East 

0.00 

45.6 

-1.39 

0.10 

Floor 

1.00 

48.3 

-1,87 

0.10 
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9.5  Results 

Figures  9.5.1  and  9.5.2  are  exactly  analogous  to  the 
coi're  spending  figures  of  Solution  V.  The  curve  for  NBS 
Test  5  is  omitted,  since  it  is  almost  identical  with  that 
in  Figure  8.5.3. 

These  results  were  also  obtained  with  the  aid  of  the 
digital  computer,  at  approximately  the  same  running  time. 

The  values  of  the  curve-fitting  constants  used  in  the 
'jomputations  are  listed  in  Table  9.5.0. 


9.6  Discussion 

Inspection  of  Figures  9.5.1  and  9.5.2  shows  the  same 
general  results  concerning  thermal  properties  found  in 
Solution  V,  except  that  Curves  2  and  3  of  Figure  9.5.2  do 
cross  near  the  end  of  the  test,  indicating  that  the  earth 
has  begun  to  take  over  a  little  earlier  than  in  Solution 
V.  Since  Solutions  V  and  VI  were  derived  on  assumptions 
identical  except  for  the  initial  temperature  function  u 
wide  \Tarlations  in  results  were  not  expected. 

Solution  VI  represents  the  last  investigation  of 
one-dimensional,  semi-infinite  region,  constant  heat  flux 
shelter  models.  A  general  discussion  and  comparison  of 
all  the  solutions  presented  in  this  report  follows  in  the 
next  section. 
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10.  GENERAL  DISCUSSION 

The  purpose  of  this  section  is  to  discuss  the  aignlfi 
cance  of  the  assumptions  made  in  the  individual  solutions 
for  four  basic  parameters,  namely: 

1)  geometry 

2)  composition  of  region 

3)  initial  temperature 

4)  thermal  properties, 
and  also  the  questions  of 

5)  shelter  irregularities 

6)  treatment  of  the  complete  problem. 


10.1  Geometry 

The  first  three  solutions  provide  a  means  of  comparing 
the  effects  of  the  basic  geometries  considered,  since  they 
differ  only  In  this  respect.  Figures  10,1,1,  .2,  and  .3 
show  surface  temperatures  versus  time  predicted  by  Solu¬ 
tions  I,  II,  and  III  for  NBS  Tests  3,  4,  and  based  on 
the  preferred  earth  thermal  properties.  The  equivalent- 
radii  used  in  Solutions  II  and  III  were  determined  on  the 
basis  of  equal  surface  areas.  Since  Solution  X,  unlike 
the  other  two,  represents  a  one-dimensional  model  and  is 
applicable  only  to  the  individual  surfaces,  it  could  still 
be  included  in  this  comparison  by  assuming  all  six  surfaces 
to  be  in  one  plane  and  using  the  same  overall  average 
values  of  heat  flux  and  initial  temperature. 

It  can  be  seen  from  these  figures  that  Solution  III, 
based  on  the  cylindrical  approximation,  produces  the  roost 
accurate  predictions  for  two-week  periods.  As  stated  in 
Section  6.3*  this  solution  might  still  be  best  for  three- 
week  periods  and  even  longer  with  increasing  size.  In  the 
case  of  large  shelters,  because  of  their  shape  and  shallow 
depth,  spherical  approximations  would  not  be  applicable. 
With  increasing  size,  the  one-dimensional  solutions  become 
Increasingly  better. 
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10.2  Composite  Region 

In  order  to  examine  the  significance  of  the  composite 
assumption,  it  is  necessary  to  compare  Solutions  I  and  IV. 
This  comparison,  however,  is  complicated  by  the  question 
of  which  set  of  thermal  properties  in  Solution  I  corres¬ 
ponds  to  which  two  sets  in  Solution  IV.  Figure  10.2.1  pro¬ 
vides  a  comparison  on  the  basis  of  the  preferred  concrete 
and  earth  thermal  properties.  The  west  wall  region  of  KBS 
Test  3  was  chosen  as  a  typical  case.  This  figure  is  a 
reproduction  of  Curves  2,  3,  and  4  of  Figure  4.4.2,  and 
Curve  1  of  Figure  7.5.2. 

Inspection  of  Figure  10.2.1  shows  Solution  IV,  Curve 
1,  lying  between  two  cases  of  Solution  I,  Curves  2  and  4. 
The  composite  assumption  lowers  the  predicted  temperature 
toward  Curve  4  for  small  values  of  time,  when  concrete 
dominates,  and  raises  the  prediction  toward  Curve  2  for 
large  values  of  time,  when  the  earth  properties  are  of  more 
significance.  This  effect  is  more  clearly  observed  by  com¬ 
paring  Solution  IV  with  Solution  I,  Curve  3,  which  is  based 
on  an  Intermediate  set  of  thermal  properties. 

It  can  be  seen  that  some  Intermediate  combination  of 
thermal  conductivity  and  diffusivity  will  translate  Curve 
3  into  a  position  of  optimum  agreement  with  Curve  1,  but 
never  into  complete  agreement,  since  the  times  under  con¬ 
sideration  are  neither  small  enough  to  approximate  the 
whole  region  by  concrete  nor  large  enough  to  approximate 
it  by  earth.  This  result  may  or  may  not  be  true  for  the 
cylindrical  and  spherical  models.  While  it  seems  likely 
that  the  effect  of  the  concrete  would  be  less  in  these 
models,  because  of  the  increasing  mass  of  earth  per  unit 
surface  area  within  a  given  distance  from  the  surface,  any 
differences  could  only  be  determined  by  solving  the  appro¬ 
priate  composite  models. 

It  was  first  pointed  out  in  Section  7.6  that  In  the 
one-dimensional  geometry,  the  earth  does  not  dominate  the 
heat  flow  sufficiently,  after  two  weeks,  for  the  wnole 
region  to  be  treated  as  earth.  This  fact  can  be  made 
plausible  by  calculating  the  relative  amounts  of  heat 
stored  in  the  concrete  and  earth,  according  to  Solution 
IV.  The  heat  stored  in  each  medium  is  simply  the  integral 
of  the  temperature  rise  over  the  medium.  Thus,  the  heat 
stored  in  the  concrete  is: 
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and  the  heat  stored  in  the  earth 


Qa<t) 
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The  temperature  rise,  0,  and  0S ,  can  be  obtained  from  Eqs. 
7.2(11)  and  (12),  giving 
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Since  F0t  is  the  total  heat  per  unit  area  lost  at  time  t, 
it  follows  that  the  fraction  stored  in  the  concrete  is  the 
coefficient  of  F0t  in  Eq.  (3),  giving  24  percent  when 
evaluated  at  t  =  336  hours.  It  may  be  noted  that  this 
agrees  fairly  well  with  an  estimate  made  in  reference  [1] , 
based  on  measured  temperatures  in  NBS  Test  3-  That  estimate 
established  the  fraction  of  heat  stored  in  all  the  concrete, 
including  the  shielding  wall,  at  about  26  percent  of  the 
total  heat  stored  within  4  feet  of  the  shelter. 
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10.3  Initial  Temperature 

Figure  10.3.1  shows  surface  temperatures  versus  time 
predicted  by  Solutions  IV,  V,  and  VI  for  the  north  wall 
region  of  NBS  Test  3>  based  on  the  preferred  sets  of  thermal 
properties.  Figure  10.3.2  shows  similar  curves  for  the  west 
wall  region.  These  figures  facilitate  a  comparison  between 
observed  performance  and  three  one-dimensional  solutions 
which  differ  only  in  the  form  of  the  Initial  temperature. 

Figures  10.3.1  and  10.3.2  show  the  effect  of  the  par¬ 
ticular  temperature  distribution  which  existed  at  the  start 
of  Test  3  on  the  predicted  surface  temperatures.  The  assumed 
temperature  at  the  shelter  surface  for  Solutions  V  and  VI 
were  as  much  as  5.6°F  lower  than  the  average  value  used  in 
Solution  IV,  resulting  in  generally  lower  predictions  from 
these  latter  solutions  for  the  first  6  to  8  days  of  the  test. 
However,  the  predicted  surface  temperature  at  the  end  of 
14  days  was  higher  from  Solution  VI  than  from  Solution  IV 
and  still  higher  from  Solution  V. 

Solution  V  was  derived  with  the  simplest  exponential 
approximation  to  the  initial  temperature  on  the  assumption 
that  the  region  beyond  4  feet  would  have  little  effect  on 
the  surface  temperature  up  to  14  days.  The  magnitude  of 
the  difference  in  the  predicted  temperatures  at  14  days’ 
time  from  Solutions  IV  and  V,  however,  suggested  that  the 
effect  of  the  earth  temperature  beyond  4  feet  might  be  sig¬ 
nificant.  It  is  probable  that  the  initial  temperature  more 
than  4  feet  from  the  shelter  surface  for  Test  3  decreased 
as  indicated  in  Figure  9.4.1,  [l].  The  effect  of  this  latter 
assumption  regarding  the  initial  earth  temperature  beyond 
4  feet  is  seen,  from  a  comparison  of  Solutions  V  and  VI,  to 
be  a  decrease  in  shelter  surface  temperature  of  a  little 
more  than  0.5°F  after  14  days. 

Solutions  V  and  VI  represent  successive  refinements 
in  the  assumptions  regarding  the  initial  earth  temperature. 

It  appears  that  the  effect  of  these  refinements  is  not  large 
after  the  first  week  of  occupancy.  Since  a  more  nonuniform 
initial  temperature  distribution  than  that  observed  for  Test 
3  is  not  likely  to  occur  in  undisturbed  earth,  the  Initial 
temperature  can  be  considered  uniform  for  design  purposes, 
unless  the  first  week  of  occupancy  is  of  comparable  import¬ 
ance  to  the  second  week. 
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Since  the  composite  nature  of  the  region  and  the  non- 
uniform  initial  temperature  have  been  taken  into  account 
in  Solutions  V  and  VI,  the  steepness  of  the  curves  in 
Figures  10.3.1  and  10.3.2  after  the  first  week  is  due 
chiefly  to  the  one-dimensional  geometry. 

10.4  Thermal  Properties 

It  was  observed  in  Sections  4.6  and  7.6,  in  discussing 
one-dimensional  Solutions  I  and  IV,  that  higher  thermal 
properties  for  the  conduction  region  produced  better  agree¬ 
ment  between  predicted  and  observed  temperatures  in  most 
cases.  The  north  wall  region  was  an  exception  to  this  trend. 
Since  heat  flows  in  one  direction  only  in  a  one-dimensional 
model,  artificially  high  thermal  properties  were  required 
to  produce  reasonable  temperature  predictions. 

It  seems  probable  that  the  NBS  family  shelter,  without 
widely  different  width,  length  and  height,  would  be  best 
represented  by  a  spherical  cavity,  especially  If  it  were 
located  deep  in  the  earth.  However,  the  comparison  of 
observed  and  predicted  surface  temperatures  by  Solution  II, 
(See.  Fig.  5.4.1  and  5.4.2),  indicate  that  the  heat  conduc¬ 
tion  region  around  the  actual  shelter  was  a  less  effective 
heat  sink  than  that  assumed  for  the  spherical  model  during 
both  the  summer  and  winter  tests.  This  was  probably  due  to 
the  effect  of  the  relatively  thin  layer  of  earth  over  the 
shelter  and  the  exposure  to  atmospheric  conditions  at  the 
earth's  surface  and  also  to  the  difference  in  shape  of  the 
cavities.  Thus,  artificially  low  thermal  properties  were 
required  for  Solution  II  to  best  approximate  the  observed 
shelter  surface  temperatures. 

Since  the  cylindrical  model  provided  a  larger  heat  sink 
per  unit  shelter  surface  area  than  the  one-dimensional  model, 
and  a  smaller  heat  3ink  per  unit  shelter  surface  area  than 
the  spherical  model,  it  tended  to  produce  the  best  agreement 
between  predicted  and  observed  surface  temperature  with  the 
preferred,  or  most  probable,  thermal  properties  of  the  heat 
conduction  medium. 

In  NBS  Tests  3  and  4,  significant  absorption  by  the 
concrete  walls  of  existing  condensation  would  have  raised 
the  thermal  properties;  however,  this  possibility  was  mini¬ 
mized  by  the  application  of  two  coats  of  rubber-base  paint. 

In  NBS  Test  5,  due  to  the  melted  snow,  the  thermal  properties 
of  the  earth  could  have  been  higher  than  the  preferred  set, 
but  this  was  not  indicated  by  the  diffusivity  determination 
discussed  in  Section  3.3. 
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The  models  used  for  this  analysis  implicitly  assumed 
that  the  heat  conduction  in  the  concrete  and  earth  was  linear 
in  nature,  i.e.,  the  thermal  properties  did  not  vary  with 
temperature,  so  the  classical  linear  heat  flow  theory  was 
applicable.  While  it  may  be  that  the  temperature  dependence 
of  the  theraial  properties  could  be  neglected  for  the  gradi¬ 
ents  encountered  in  these  tests,  these  gradients  nevertheless 
may  have  caused  some  moisture  migration  in  the  soil.  This 
migration  would  have  acted  both  to  transfer  sensible  and 
latent  heat  directly,  and  to  change  the  thermal  properties, 
which  are  much  more  sensitive  to  moisture  changes  than  to 
temperature  changes  £7] .  The  result  is  a  problem  of  com¬ 
bined  heat  and  moisture  flow  which  involves  coupled,  non¬ 
linear,  partial  differential  equations,  amenable  to  solution 
by  numerical  analysis  [13J  . 

Although  a  few  measurements  of  moisture  content  of  the 
earth  were  made  during  the  course  of  the  five  NBS  tests, 
the  data  did  not  indicate  whether  significant  moisture 
transfer  occurred.  It  probably  was  not  sufficient  to  cause 
significant  changes  in  the  thermal  properties  of  the  earth. 
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10.5  Shelter  Irregularities 

Throughout  this  analysis ,  both  heat  conduction  media 
were  assumed  to  be  homogeneous,  and  in  the  case  of  the 
earth.  Infinite  in  extent.  These  conditions  were  not  fully 
met  by  the  NBS  test  shelter  and  site,  due  to  certain 
irregularities  in  the  walls  and  earth.  The  deviations  for 
each  region  are  listed  as  foDlows: 

1)  north  wall  region  —  a  concrete  road  was  located 
about  ten  feet  from  the  wall,  which  might 
actually  have  acted  as  a  heat  source,  storing 
heat  from  the  sun.  This  additional  heat  source 
would  tend  to  make  the  surface  temperature  rise 
faster,  requiring  lower  thermal  properties  to 
make  predictions  agree  with  observations  (cf. 
Sections  **.6  and  7.6).  The  hatchway  e*it, 
adjacent  to  this  wall,  might  have  affected  the 
distribution  of  heat  flux  over  the  wall.  The 
shielding  wall  might  also  have  affected  the  heat 
flow  distribution,  by  interfering  with  the  radia¬ 
tion  component  of  the  internal  heat  exchange. 

This  wall  also  contained  the  metal  exhaust  pipe, 
which  would  tend  to  conduct  some  heat  directly 

to  the  air.  It  is  for  these  reasons  that  this 
region  was  designated  the  least  typical. 

2)  west  wall  region  _  the  south  half  of  this  wall  was 
probably  exposed  to  more  heat  than  the  ncx-th  half, 
due  to  the  arrangement  of  simulated  occupants  and 
shielding  wall,  but  since  it  faced  uphill  on  a  site 
of  12  percent  grade,  it  was  exposed  to  more  earth 
than  any  other  wall.  This  region  was  designated 
the  most  typical. 

3)  south  wall  region  —  the  ventilating  air  supply 
system  used  during  the  NBS  tests  was  located  in 
the  shielding  wall  and  aimed  at  the  south  wall. 
Furthermore,  three  simulated  occupants  were  located 
near  this  wall.  This  arrangement  possibly  affected 
the  heat  flux  distribution.  A  spare  ventilating 
air  supply  and  exhaust  system,  installed  in  this 
wall,  would  tend  again  to  conduct  some  heat  directly 
to  the  air. 
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4)  east  wall  region  —  this  wall  faced  downhill  and 
thus  was  exposed  to  less  earth  than  any  other  wall. 
In  fact,  the  original  grade  line  was  about  one 
foot  below  the  top  edge  of  this  wall.  It  was 
consequently  the  one  most  affected  by  weather  con¬ 
ditions,  Also,  the  hatchway  portion  of  this  wall 
received  less  heat  than  the  rest. 

5)  floor  region  —  this  region  had  a  truly  infinite 
exposure  of  earth.  However,  the  hatchway  portion 
of  the  floor  received  less  heat,  and  the  arrange¬ 
ment  of  the  simulated  occupants  and  the  supply 
air  stream  might  have  affected  the  heat  flux 
distribution. 


10.6  Shelter  Environment 

In  this  report,  the  conduction  problem  has  been  treated 
as  separate  from  the  rest  of  the  shelter-environment  system, 
using  the  known  heat  flux  at  the  surfaces.  While  this  pro¬ 
cedure  is  not  invalid,  it  is  incomplete  in  the  sense  of 
omitting  that  part  of  the  heat  transfer  process  between  the 
internal  heat  sources  and  the  walls.  A  more  complete  treat¬ 
ment  would  take  Into  account'Hhese  sources,  namely,  venti¬ 
lation  and  occupants.  In  such  a  problem,  it  would  not  be 
possible  to  treat  each  surface  separately,  as  was  done  In 
the  one-dimensional  solutions,  because  the  heat  losses  to 
all  the  surfaces  are  interrelated  through  the  common  sources. 
In  that  case,  uhe  simplest  approach  would  be  to  use  a  single 
symmetrical  model  such  as  the  cylinder  or  spnere,  and 
according  to  Section  10.1,  the  cylinder  is  the  better  choice. 

A  cylindrical  solution  with  heat  source  boundary  condi¬ 
tions  is  available  [  143,  and  has  been  expanded  and  applied 
to  the  NBS  shelter.  However,  the  composite  problem  has  not 
been  found  in  a  literature  search,  and  the  composite  assump¬ 
tion  may  be  important  as  seen  in  Section  10.2. 
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When  the  conduction  problem  is  treated  separately, 
the  roof  region  can  be  accurately  represented  by  a  one- 
dimensional  finite  layer,  with  boundary  conditions  at  both 
ceiling  and  ground  surfaces.  (Three  such  solutions  have 
been  developed,  but  are  not  included  in  this  report.) 
Because  of  the  proximity  of  the  ground  surface,  the  heat 
flow  from  the  four  walls  is  chiefly  two-dimensional,  ex¬ 
cept  near  the  floor.  The  floor  is  exposed  to  the  most 
earth  of  all  the  surfaces,  and  heat  flow  in  that  region 
is  truly  three-dimensional.  These  facts  suggest  that  the 
overall  effect  might  be  two-dimensional,  adding  weight  to 
the  choice  of  cylindrical  geometry  as  the  best  simple 
model. 

It  should  be  noted  that  none  of  the  one-dimensional 
solutions  in  this  report  are  really  applicable  to  the 
floor  region,  because  of  its  three-dimensional  behavior. 
However,  3ince  the  floor  region  was  approximately  semi¬ 
infinite,  the  predicted  temperatures  were  included  in  the 
tables  for  the  sake  of  completeness. 

A  variety  of  more  elaborate  models  might  be  considered 
which  would  explicitly  take  into  account  the  asymmetry 
caused  by  the  ground  surface.  An  example  might  be  a  closed 
hemispherical  cavity  convex  downward  located  at  a  finite 
depth  within  a  semi-infinite  region.  This  model  would 
probably  exhibit  a  more  desirable  heat  flow  pattern,  approx' 
imately  one-dimensional  in  the  roof  region,  two-dimensional 
at  the  sides,  and  three-dimensional  at  the  bottom.  Models 
of  this  type  are,  of  course,  mathematically  complex. 
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11,  CONCLUSIONS 

On  the  basis  of  the  analyses  covered  by  this  report* 

It  is  concluded  that  a  cylindrical  model  provides  the  best 
predictions  of  the  surface  temperatures  in  a  small  under¬ 
ground  shelter  with  limited  earth  cover*  probably  because 
it  most  nearly  compensates  for  the  effect  of  the  finite 
region  over  the  cavity.  It  is  probable  that  a  spherical 
model  would  be  better  if  the  small  shelter  were  considerably 
deeper  in  the  earth.  The  one-dimensional  model  was  less 
effective  in  predicting  temperatures  in  such  a  small  shelter, 
but  would  probably  be  an  adequate  representation  for  larger 
shelters  in  which  comer  effects  were  relatively  of  less 
importance. 

The  composite  nature  of  a  concrete  shelter  surrounded 
by  earth  was  shown  to  have  a  small  effect  on  shelter  surface 
temperatures.  The  higher  thermal  properties  of  concrete  a3 
compared  to  earth  caused  the  shelter  surface  temperatures 
to  be  a  little  lower  than  if  the  entire  conduction  region 
was  comprised  of  earth.  This  effect  was  slightly  greater 
during  the  first  week  of  occupancy  than  during  the  second. 

A  heat  conduction  region  comprised  entirely  of  concrete 
might  reduce  surface  temperatures  as  much  as  3°P  in  a  small 
shelter  as  compared  to  a  shelter  with  only  moist  earth  as  a 
heat  sink.  This  disparity  would  be  even  greater  if  the  earth 
medium  was  essentially  dry. 

The  initial  distribution  of  temperature  from  the  shelter 
surface  to  a  distance  of  4  feet  into  the  surrounding  medium 
was  found  to  make  an  appreciable  difference  In  the  surface 
temperature  during  the  first  week  of  occupancy,  but  the  ef¬ 
fect  was  of  little  significance  after  two  weeks. 

The  solutions  described  in  this  report  do  not  take  into 
account  the  effect  of  moisture  migration  on  heat  conduction, 
but  these  effects  were  believed  to  be  small  in  the  test  snel- 
ter. 


Further  efforts  toward  fulfilling  the  overall  objectives 
of  this  analytical  study  should  begin  with  models  that  take 
into  acco  'nt  heat  generation  and  ventilation  of  the  shelter, 
the  finite  roof  cover,  and  the  composite  nature  of  the  shel¬ 
ter  surroundings. 
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APPENDIX  A 

INITIAL  TEMPERATURE  DISTRIBUTION 
IN 

SHELTER  CONCRETE 
AND 

SURROUNDING  EARTH  MEDIUM 
FOR 

FIVE  NBS  SHELTER  TESTS 
CONDUCTED  BETWEEN 
August  13,  1959  and  April  8,  i960 
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NBS  TEST  3 

INITIAL  TEMPERATURE  DISTRIBUTION  IN 
HEAT  CONDUCTION  REGION  ON  SEPT.  13.1959 


N8S  TEST'S 

INITIAL  TEMPERATURE  DISTRIBUTION  IN 
Lier  at  AAkim  irTIftkl  PCfllAW  DM  25.  I960 
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APPENDIX  B 
Heat  Flux 

The  mathematical  solutions  in  this  report  are  based 
on  a  constant  ( time-independent)  heat  flux,  applied  uni¬ 
formly  over  the  entire  bounded  surfaces  of  each  mathe¬ 
matical  model.  These  are  conditions  which  would  appear 
more  readily  applicable  to  a  sealed-up  shelter,  and  which 
would  fail  to  describe  the  case  of  ventilated  shelters, 
whose  heat  flux  would  not  be  expected  to  be  either  con¬ 
stant  or  uniformly  applied.  Thus,  in  order  to  deal 
properly  with  the  constant  flux  parameter  appearing  in 
each  solution,  it  was  necessary  to  extract  numerical 
values  from  the  experimental  data  based  on  methods  con¬ 
sidered  to  be  compatible  with  the  constant  flux  assumption. 
Three  methods  were  considered  and  each  is  discussed  below, 
one  of  which  required  the  derivation  of  a  new  solution. 

Method  1  Constant  Heat  Flux 

Each  solution  based  on  this  assumption  required  that 
a  single  value  of  heat  flux  be  used  for  all  values  of 
time.  This  value  could  reasonably  be  selected  as  the 
overall  average  value,  shown  in  Table  B.l  at  336  hours, 
and  which  is  close  to  the  experimental  value  observed 
during  approximately  the  seventh  day  of  each  test,  for 
any  surface.  However,  the  observed  performance  diverges 
from  this  value  for  other  times. 

Method  2  Cumulative  Average  Heat  Flux 

This  is  the  method  used  throughout  the  report.  The 
values  of  heat  flux  were  cumulative  averages,  values  which 
are  essentially  time-dependent,  averaged  over  the  interval 
from  o  to  t,  and  considered  to  be  constant  in  that 
interval.  The  smaller  intervals  gave  a  better  overall 
approximation  than  Method  1  because  the  values  thus  ob¬ 
tained  were  distributed  such  as  to  incorporate  the  ob¬ 
served  exponential  decay  feature.  Even  though  this 
method  was  obviously  a  deviation  from  the  models,  it  was 
nevertheless  considered  to  be  a  reasonable  compromise 


T 

i,. 

i* 
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between  the  mathematical  assumption  and  the  experimental 
observations.  The  predicted  curve  determined  on  this  basis 
is  then  actually  a  curve  drawn  through  points  obtained  by 
solving  each  solution  eight  separate  times,  each  time  vary¬ 
ing  the  heat  flux  in  accordance  with  the  cumulative  average 
values  listed  in  Table  B.l. 

Method  3  Prescribed  Heat  Flux 

This  method  would  have  required  the  derivation  of 
each  solution  on  the  basis  of  an  actual  time-dependent 
heat  flux  boundary  condition,  but,  as  previously  mentioned, 
this  was  considered  to  be  beyond  the  scope  of  this  report. 
However,  in  order  to  determine  the  effects  of  the  mathe¬ 
matical  simplifications  involved  in  Method  1,  '*nd  of  the 
"compromise"  made  in  Hethod  2,  Solution  I  was  reworked  by 
replacing  Eq.  4.2(3)*  the  constant  flux  condition,  by 


-K  f§  =  a  +  be-mt  a: 

a  prescribed  flux  condition  approximating  the  experi¬ 
mental  situation  of  exponentially  decaying  heat  flux. 
This  solution  is  Identified  as  Solution  IA. 

The  temperature  rise  0  at  the  interior  surface  x 
is  then 


°  =  ¥ Iff  ♦ ! te? 


where  D(^/mt)  represents  Dawson’s  integral,  equal  in  this 
case  to  the  imaginary  part  of  the  complex  error  function. 
Both  functions  are  tabulated  in  [15] • 
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Figures  B.l  and  B.2  illustrate  predicted  surface 
temperatures  obtained  from  Methods  1*  2,  and  3»  together 
with  the  observed  data,  for  che  north  and  west  wall  regions 
of  NBS  Test  3«  In  both  cases,  the  prescribed  flux  solution. 
Solution  IA,  yielded  results  in  closer  agreement  with  the 
observed  temperatures  at  fourteen  days.  Method  1  is  gen- 
erally  better  for  the  first  seven  days  because  the  constant 
flux  value  used  was  lower  than  or  about  equal  to  the  ob¬ 
served  values  for  those  times.  At  the  end  of  fourteen 
days,  the  deviations  in  temperature  between  Methods  2  and 
3  were  less  than  1°P  for  the  west  wall  and  practically 
negligible  for  the  north  wall.  For  times  larger  than  336 
hoursj  results  from  Methods  2  and  3  will  tend  to  merge, 
that  trend  being  already  visible  in  Figure  B.l,  so  that 
the  greatest  deviation  is  probably  encountered  at  times  in 
the  neighborhood  of  fourteen  days  or  less. 

A  similar  improvement  using  Bq.(l)  could  be  made  to 
the  other  solutions  in  this  report,  but  the  likelihood  of 
uncovering  deviations  much  wider  than  those  illustrated 
by  this  solution  appears  doubtful. 
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TABLE  B.l 


CUMULATIVE  AVERAGE  VALUES  OP 
HEAT  FLUX  F0 

Test  3 


Time 

North 

West 

South 

East 

Roof 

Floor 

2  4 

2.06 

4.41 

8.32 

4.78 

3.60 

5.24 

48 

2.16 

4.02 

7.09 

4.44 

3.34 

5.17 

72 

2.15 

3.70 

6.21 

4.13 

3.18 

5.12 

96 

2.02 

3.56 

5.57 

3.90 

2.97 

5.02 

120 

1.93 

3.41 

5.10 

3-71 

2.86 

4.90 

168 

1.85 

3.25 

4.50 

3.46 

2.80 

4.66 

240 

1.76 

3.04 

3.99 

3.22 

2.70 

4.50 

336 

1.65 

2.84 

3-53 

3.01 

2.36 

4.24 

Test  4 


Time 

North 

West 

South 

East 

Roof 

Floor 

24 

2.65 

4.78 

5.49 

4.71 

3.96 

4.70 

48 

2.63 

4.61 

4.97 

4.42 

3.45 

4.51 

72 

2.53 

4.42 

4.64 

4.18 

3.15 

4.40 

96 

2.48 

4.34 

4.40 

4.02 

2.93 

4.35 

120 

2.45 

4.28 

4.26 

3.90 

2.77 

4.32 

168 

2.31 

4.10 

3.98 

3.74 

2.41 

4.23 

240 

2.07 

3.87 

3.63 

3.42 

2.27 

4.0^ 

336 

1.83 

3-69 

3.35 

3.18 

2.36 

3.90 

Test  5 


Time 

North 

West 

South 

East 

Roof 

Floor 

24 

1.81 

4.19 

6.38 

5.74 

4.57 

1.73 

48 

1.82 

3.99 

6.08 

5.35 

4.27 

1.74 

72 

1.85 

3.93 

5.96 

5.26 

4.16 

1.76 

96 

1.89 

3.89 

5.88 

5.18 

4.01 

1.78 

120 

1.90 

3.84 

5.79 

5.10 

3.80 

1.81 

168 

2.02 

3.91 

5.80 

5.09 

3.44 

1.95 

240 

2.09 

3.90 

5.74 

4.97 

2.99 

2.08 

336 

2.06 

3.86 

5.65 

4.76 

2.80 

2.19 

SOLUTION  IA  -  NBS  TEST  3 
PREDICTED  SURFACE  TEMPERATURES 
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APPENDIX  C 
Numerical  Example 

The  example  worked  out  below  for  Solution  IV  is  the 
type  of  computation  necessary  to  obtain  numerical  results 
from  any  of  the  first  four  solutions  contained  in  this 
report . 

To  compute  the  Interior  surface  temperature  Uj ,  at  x  =  -l, 
for  the  north  wall  of  NBS  Test  3  at  336  hours,  Eq.  7*2(13) 
is  manipulated  as  follows. 

Rewriting  the  summation  portion  of  the  equation  in  a 
slightly  more  convenient  form  yields 


2<-Dn  @n  { 


2  ierfc  -2 


-20  ierfc 


n+1 

T8 


} 


(1) 


The  factor  of  2  was  brought  Inside  the  summation  sign  because 
the  integral  error  function  of  the  first  order  is  usually 
tabulated  in  this  form.  Numerical  values  for  this  function, 
as  well  as  for  those  of  higher  order  appearing  in  Solution  III, 
are  tabulated  In  [15]  and  in  Appendix  II  of  [10]. 

Expanding  summation  (l)  to  produce  sufficient  convergence 
for  the  desired  accuracy  for  the  particular  values  of  time, 
wall  thickness,  and  thermal  properties  involved  will  require, 
in  this  example,  the  first  four  verms,  and  a  convenient 
arrangement  for  listing  these  terms  is  shown  below.  At 
n=o,  1,  2,  3***<  the  summation  reduces  to  the  following: 


n  =  o,  2  ierfc  (o)  -  20  ierfc  ["" 


-20  ierfc  \  -i  1  +  20*ierfc  f 

L  J  L  Tf 


n  =  1, 


* 
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T 


n  =  2,  2ga 


ierfc  f  -4  1  -  2@s  ierfc  T  *4  1 

>-  -J  L 


n  *  3 ,  -2(3*  ierfc  [  ~r  J  +  negligible  term 


Eq.  7.2(13)  can  be  rewritten,  combining  the  above 
i/erms,  to  read 


The  following  substitutions  can  now  be  made  in  the 
expression  to  obtain  the  required  numerical  result; 


Prom  Table 

u0 

- 

71.7 

Prom  Table  B.l, 

Fo 

■s 

1.63 

From  Table  7-5.1> 

K, 

— 

1.15 

ai 

— 

0.036 

a 

= 

0.767 

Thickness  of  concrete. 

l 

s 

0.667 

above 


Therefore, 


T 


04 1 

~7 


=  27.216 


8  =  =  -0.1316 

C~!-l 


and 
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Substituting  In  (2),  we  get 


ux  =  71.7  +  (5.2169)  {  1.1284  +  0.2632(0.7862) 

1.15 


+  0.0346(  0.5235)  +  0.0046(  0.3321)  4 


} 


uj  =  71.7  +  6.76  =  78.46  or  78.5°F, 
the  predicted  value  shown  in  Table  7*5.1. 

It  should  be  noted  that  Solutions  V  and  VI  are  too 
cumbersome  for  the  type  of  manipulation  shown  in  the  above 
example,  especially  when  many  computations  are  required. 
Even  when  convergence  is  fairly  rapid,  the  use  of  a  digital 
computer  is  almost  a  necessity.  As  indicated  in  Sections 
8-5  and  9.5,  these  two  solutions  were  programmed  in 
Fortran  language  for  use  on  the  IBM  7090  computer,  and 
these  programs  are  available  from  the  Mechanical  Systems 
Section  of  NBS. 
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